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Abstract 

We present an algorithm to find an Euclidean Shortest Path from a source vertex s to a sink 
vertex t in the presence of obstacles in 5ft 2 . Our algorithm takes 0(T + m(lgm)(lgn)) time and 
0(n) space. Here, 0{T) is the time to triangulate the polygonal region, m is the number of 
obstacles, and n is the number of vertices. This bound is close to the known lower bound of 
0(n + mlgm) time and 0(n) space. Our approach involve progressing shortest path wavefront 
as in continuous Dijkstra-type method, and confining its expansion to regions of interest. 

1 Introduction 

The shortest path problem in 3fJ d is that of finding a shortest route from one point to another among 
the presence of obstacles. Even in K 3 under the Euclidean metric, it is not even known whether 
the shortest path problem in the presence of polyhedral obstacles is in NP though the problem 
has been shown to be NP-hard. This paper considers the case in K 2 . The Euclidean shortest path 
problem in a polygonal region is one of the oldest and best-known in computational geometry due 
to its various applications. Mitchell [13J provides an extensive survey of research accomplished in 
determining shortest paths in polygonal and polyhedral domains. 

We assume that the domain is defined by a simple polygon having m obstacles comprising a 
total of n vertices. There are two fundamentally different approaches in solving this problem: the 
visibility graph method, and the wavefront method. These approaches assume a triangulation of 
the domain, which can be accomplished in 0(n + m(lgm) 1+e ) using the algorithm from Bar- Yehuda 
and Chazelle pQ. 

The visibility graph method is based on constructing a graph whose nodes are the vertices of 
the obstacles and the edges are pairs of mutually visible vertices. Welzl [17] provides an algorithm 
for constructing the visibility graph with n line segments in 0(n 2 ) time. Ghosh and Mount [3], and 
Kapoor and Maheshwari [9] found an algorithm to construct the visibility graph of time complexity 
0(nlgn + E), where E is the number of edges in the graph. Applying Dijkstra-type algorithm on 
this graph, one can determine a shortest path in 0(nlgn + E). Unfortunately the visibility graph 
can have Q(n 2 ) edges in the worst case, so any shortest path algorithm that depends on an explicit 
construction of the visibility graph will have a similar worst-case running-time. 

Storer and Reif |15J presented 0(T + mn) time algorithm which constructs a data structure so 
that the shortest path from s to any point on the plane can be determined in 0(1) time. Using 
the concept of corridors, Kapoor and Maheshwari [8] presented an algorithm of time complexity 
0(m? lgn + nlgn). 
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The second approach used by Hershberger and Suri [3], Mitchell [TT], Mitchell [12], and Kapoor 
[BJ [7] gave algorithms to find a shortest path by expanding a wavefront from source s till it reaches 
the destination t. This approach seems inherently more geometric than the graph-theoretic method 
based on visibility graphs. However, this method when directly applied does not achieve the known 
lower bound of fi(n + mlogm), and resolving this is an open problem for several years. Mitchell 
|12j gave an algorithm for computing a shortest path map, an encoding of shortest paths from s 
to all points of the plane in 0(n 3 / 2+e ) time and space. More recently, Hershberger and Suri [5] 
presented O(relgn) time algorithm. 

In this paper, we combine corridors with the wavefront approach to obtain a 0(n+m(lg m)(lg n)) 
time algorithm which uses 0(n) space for computing Euclidean shortest path among obstacles in 
whereas the Problem 21 of The Open Problems Project (TOPP) intends for a solution with 
0{n + m\gm) time using 0{n) space. We assume a model of computation where real arithmetic is 
allowed, though the results apply even when finite precision arithmetic is used (ignoring the numer- 
ical complexity of the schemes). Our algorithm proceeds by first triangulating the given polygonal 
region and then identifying the useful corridors and junctions among those triangles as in Kapoor 
et al. [8]. Then we initiate a shortest path wavefront from source and progress it as in continuous 
Dijkstra-type of method; however, to reduce the number of event points, we confine the wavefront 
to progress in regions of interest. 

Section [2] gives basic definitions, properties, and the utility of various constructs that we use 
in developing the algorithm. Algorithm outline is mentioned in Section El Section [4] gives details 
of data structures and the operations on each of them. More technical details of algorithm are 
presented in Sections El El El El and [TTJl The algorithm in terms of event point types and their 
handling is described in Section El Both the analysis and proof of correctness are spread all 
through the paper, whereas the analysis required for the overall time and space complexity analysis 
is presented in Section [TTJ Section [T2] concludes with possible generalizations. 

2 Definitions and Properties 

Let vi=vq, vi, . . . , be I points, known as vertices, in the plane. The sequence of I line segments, 
known as edges, ei=eo=VQVi, e\=v\V2, ■ ■ ■ , e;_i=w;_iWo together form a closed polygonal chain, say C. 
The polygonal chain C is simple if and only if VjeiHej+i = Uj+i and Vj j-^ieiDej = <f>. Then C with 
the region bounded by it together is known as a simple polygon, say P'. Let O = {Pi, Pi, . . . , P m } 
be the set of simple polygons interior to C s.t. Vij^jPj n Pj = <j>. Then P — (J i P, is known as the 
polygonal domain or polygon with holes, say V . The set V comprising the vertices of V is of size 
n, whereas the number of obstacles \0\ is m. We denote the Euclidean shortest distance between 
two points pi and p2 in P with d(p\,p2)- We intend to find an Euclidean shortest path between 
two points in V . Among these two points, one is termed as source s and the other is sink t. We 
consider both s and t as degenerate single point obstacles. In other words, s,t 6 O. We use the 
continuous Dijkstra's approach in finding a shortest path from s to t. 

Wavefront Progression with Triangulation 

Definition 2.1 The shortest path wavefront W(d) at distance d is the locus of points at Eu- 
clidean shortest distance d from the source s. 

Initially, the wavefront is a circle centered at s with radius e, where e is a small positive 
constant. The algorithm proceeds by expanding the wavefront in "P. As the wavefront progresses, 
it may encounter various vertices and edges of V. Let the shortest path wavefront be at distance 
d' from s. A point p £ V is considered as traversed if d(p, s) < d' . Otherwise, the point p is said 
to be untraversed. An edge e is traversed, if there exists a point p € e such that p is traversed. A 
region R is traversed, if for every point p € R, p is traversed. An edge e is defined as struck if 
there exists a point p £ e such that d{p, s) = d' and for p ^ p' and p' £ e, d{p' , s) > d' . For any 
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(a) Before applying the property (b) After applying the property 

Figure 1: Non-contiguity property of shortest paths 



vertex v of V, when the wavefront strikes v, a new arc with center v may be initiated and inserted 
into the wavefront. 

Property 2.1 A wavefront segment w(v) is a circular arc with center v € V , such that each 
point on w(v) has a shortest path to source s via v. 

Property 2.2 At any stage of the algorithm, the wavefront W(d) comprises a contiguous (abutting) 
sequence of wavefront segments, w(vi), . . . , w(vk) for some k < n. 

The algorithm halts when the wavefront strikes t. Suppose w(vk) is the arc that struck t, 
w{vk-i) is the arc that struck Vk, • • • , w(s) is the arc that struck v\. Then our algorithm outputs 
the shortest path from s to t which comprises of line segments svi, i»ii?2, • • • , v^-iv^ and v^t. Adding 
the Euclidean distances along these line segments yields the shortest distance from s to t. 



Lemma 2.1 Consider shortest paths from source s to two points pi and pj in V . Suppose the 
interior of a line segment in shortest path from s to pi intersects with the interior of a line segment 
in shortest path from s to pj. Then there always exists a shortest path from source s to pi (resp. s 
to pj) so that the interior of no line segment in shortest path from s to pi (resp. s to pj) intersects 
with the interior of line segments in the given shortest path from s to pj (resp. s to pj). This 
property is termed as non-crossing property of shortest paths. 

Proof: 

Consider the point p at which the interior of line segments p[p'l and p'jP'i respectively belonging 
to shortest paths from s to pi and pj intersect. First, note that there exists at least two shortest 
paths from s to p: one via p\ and the other via p'-. Replacing the given shortest path to p'- via p\ 

with p "p and pp'j , and replacing the given shortest path to p" via p'j with p'-p and pp\ yields two 
shortest paths that do not cross with each other at p. Since this local operation reduces the number 
of interior intersection of line segments in both the shortest paths by one, repeatedly applying this 
operation at every such intersection point yields the required. □ 

To guide the wavefront progression, we triangulate the polygonal domain using the algorithm 
by Bar- Yehuda and Chazelle pQ. First we triangulate V — {s,t}, and obtain a triangulation. We 
locate s in a triangle T = (vi,V2,vs) and introduce triangulation edges sv\, SV2, sv^); similarly, after 
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Figure 2: Corridors and Junctions 



locating t, we introduce triangulation edges to yield a resultant triangulation, denoted with T. 

A sequence of wavefront segments in the wavefront, not necessarily contiguous, together are 
termed as a section of wavefront. We define a line segment I to be reachable by a section of 
wavefront, if there exists a point p on I and a point w on W such that the interior of line segment 
pw does not intersect any untraversed vertex/edge in V . Typically, we are interested in reachable 
edges of triangulation T ■ The wavefront is progressed based on its interaction with the reachable 
edges, and the interaction among the wavefront segments. 

The algorithm is event based. Let W(d) be the wavefront. Primarily, the event points can 
be categorized into two: finding d' > d so that the W(d') strikes a reachable triangulation edge; 
finding d" > d so that the W(d") at which two wavefront segments intersect. The events occur 
as the wavefront progresses and are maintained in a min-heap, with the corresponding shortest 
distance at which the event occurs (d! or d") as the key. The former causes the updates to the 
set of reachable edges, and the new segments may possibly be incorporated into the wavefront. 
The latter could cause some wavefront segments to change shape, some to be removed from the 
wavefront, and may change the set of reachable edges. 

With this approach, there are 0(n) vertices from each of which a wavefront segment could 
be initiated; these 0(n) wavefront segments could interact with 0{n) triangulation edges, causing 
quadratic time complexity in terms of number of vertices in V. 

Corridors and Junctions 

We intend to reduce the number of edges with which the wavefront interacts. This we ac- 
complish by exploiting the structure in the triangulation and obtain a coarser data structure. The 
number of elements in the new data structure are 0(m) and the wavefront may strike 0(m) entities 
rather than 0(n), hence is an improvement. 

The coarser structural elements obtained from triangulation are termed as corridors and junc- 
tions. We outline their descriptions from Kapoor et al. [10]. Consider the triangulation T and its 
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dual graph. See Fig. [2j The dual graph is divided into paths where a path is composed of maxi- 
mally connected vertices of degree two. Each such path defines a corridor formed by the sequence 
of dual regions or triangles corresponding to the vertices of the dual graph. A corridor is a region 
confined by (at most) four geometric entities - two convex chains on opposite sides, termed as cor- 
ridor convex chains; and two edges, termed as (wavefront) enter/ exit bounding edges or enter /exit 
boundaries, that are incident to both the corridor convex chains. Each of these convex chains' is a 
section of boundary of an upper or lower hull (see [2j). A junction is a triangle enclosed by at most 
three edges, each edge from a distinct corridor. A side of an edge belonging to a useful corridor 
convex chain or enter/exit bounding edge is known as a bounding edge. An end point of a bounding 
edge is known as a bounding vertex. An edge e in the dual graph is defined to be useful if there 
exists a simple path from s to t that is having a point in common with the triangles associated with 
dual edge e. Otherwise, it is useless. A corridor is defined to be useful if the triangles constituting 
the corridor contribute an useful edge in the dual graph. While partitioning V into corridors and 
junctions, we define s and t as degenerate corridors. The algorithm starts with a subdivision of the 
polygonal region into 0(m) useful corridors and junctions. 

The corridors can be classified by their structure into two types, open and closed. A corridor C 
is termed as an open corridor whenever there exists two points such that p± lies on one enter/exit 
boundary of C and P2 lies on the other enter/exit boundary edge of C, and p\ is visible to p2- See 
Fig. [3l Otherwise, a corridor is termed as a closed corridor. A closed corridor gives rise to two 
funnels each with an apex, and each funnel has two convex chains. See Fig. HI 

The advantages in progressing shortest path wavefront using corridors is two fold. First, rather 
than interacting with the 0{n) triangulation edges, it interacts with 0{m) corridor convex chains 
and corridor enter /exit bounding edges. Further, as explained latter, we exploit the coherence in 
the wavefront segments that are originated from the successive vertices along a corridor convex 
chain. 

A corridor convex chain C is reachable, if there exist a bounding edge b € C that is reachable. 
A contiguous cycle of reachable corridor convex chains and/or enter /exit boundaries is termed as 
a boundary cycle. 

As shown in Fig. El^a), ei,e2,e3 are the edges of a junction J and the corresponding 61,62,63 
are the bounding edges that are reachable from the initial wavefront segment w(s). The initial 
wavefront comprises of w(s) only. The bounding edges 61,62,63 together form the boundary cycle 
BC. The first event that occurs when the wavefront, which is a circle, strikes 61 of J. See Fig. [3](b) . 
At that stage, the bounding edge 61 in BC is replaced by convex chain C[, bounding edge 64, and 
the convex chain C'{. The resultant boundary cycle is C[, 64, C", 62, 63. As the w(s) progresses, the 
boundary cycle BC further changes as shown in Fig. EKc). In general, as the wavefront expands, 
if the just struck edge e bounds an untraversed junction J then e is replaced by the other two 
edges of J in the boundary cycle under consideration. And if the edge e is a bounding edge of an 
untraversed corridor C, then e is replaced by the other bounding edge of C and corridor convex 
chains of C. 

In the case of closed corridors, after the wavefront strikes the first apex point of the funnel, 
segments are initiated from the other apex point (provided the shortest distance to that has not 
already been determined) p, when the wavefront expands from s after a distance that equals the 
shortest distance from s to p. See Fig. HI 

As the wavefront progresses, a boundary cycle may split. See Fig. [5j This is possible when an 
edge e in a boundary cycle BC is reachable from the wavefront from both of its sides. In other words, 
e appears twice in BC as a bounding edge, and the boundary cycle splits into two when e is struck 
from either side. Consider a boundary cycle BC = 61, 62 ... 6j = (ut,Vi) . . . bj = (vi,Ui) . . . 6&, bi. 
Note that both 6j and bj represent the same edge, say e. When e is struck, the boundary cycle 
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(a) Boundary cycle bi, bi, bz (b) Boundary cycle C[, 64, C", &2, 63 




(c) Boundary cycle C(, 64, C" , b 2 ,C 3 ,b 6 , C'i 



Figure 3: Boundary cycle as the initial wavefront expands 
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Figure 4: Wavefront progression in closed corridors 



BC splits into two boundary cycles: BC{ = bi, . . . , ■ ■ ■ ,bk and BC, = . . . , 61-1 . For 

any two boundary cycles BC',BC", and bounding edges 6' <G BC',b" <G BC", if b' and 6 do not 
correspond to the same edge, then we say that BC', BC" are disjoint boundary cycles. Therefore, 
BCi and BCj are disjoint boundary cycles. The corridor bounding edges in all boundary cycles 
together are denoted with dB. A sequence of contiguous bounding edges in dB is termed as a 
section of boundary. 

Associations and I-curves 

Since there are 0{m) corridors, 0(m) junctions, and 0(n) wavefront segments, the interactions 
between the wavefront and the junction/corridor boundaries can be 0{nm). This has been per- 
ceived to be the bottleneck in implementing the wavefront method. To this end, we keep track 
of corridor convex chains, corridor enter/exit boundaries that are reachable from sections of wave- 
front, with the wavefront progress. 

Consider a corridor convex chain or enter/exit boundary g in a boundary cycle. Let W g be the 
set of wavefront segments from each of which g is reachable. Every wavefront segment w € W g is 
associated with g (or, g is associated with segment w) if and only if a point on g has shortest 
Euclidean distance to s via the center of w. The association is defined by the relation: A C Q x S, 
where Q is the set of corridor convex chains or enter/exit boundaries in the polygonal domain, and 
S is the set of waveform segments formed during the course of algorithm. Note that the wavefront 
segments in W g need not be contiguous in wavefront. 

We maintain these associations for determining the wavefront progression at which the wave- 
front strikes dB. For every section of wavefront W, given its association with a section of boundary 
BS in a boundary cycle, we need only to compute the shortest distance between BS and W, which 
is more efficient than computing the shortest distance between W and whole of dB. Again, we 
update these associations locally whenever an event changes either BS or W. The following defi- 
nition of I-curves helps in initiating and updating associations of wavefront segments or sections of 
wavefronts. 

The Voronoi diagram of a given set of points S is the partition of the plane into regions so 
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(a) Junction J before boundary split 




boundary cycle (BC , ) boundary cycle (BC 2 ) 
(b) Junction J after boundary split 

Figure 5: Boundary split 
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Figure 6: Contiguity property for segments 



that all the points interior to a region are closer to one and only one point of S. To determine the 
associations for a wavefront segment, essentially we require Voronoi regions corresponding to each 
wavefront segment. Consider any two adjacent wavefront segments, w(v a ),w(vb), in the wavefront. 
The curve that separates the Voronoi regions of v a and Vf, is termed as an I-curve(w(v a ),w(vb)). 
Every point p on I-cuive(w(v a ),w(vb)) has (at least) two shortest distance paths to s: one is via 
v a (the line segment from p to v a together with a shortest path from v a to s) and the other is via 
Vf, (the line segment from p to Vb together with a shortest path from Vb to s). In other words, I- 
cwcve(w(v a ),w(vb)) is a curve separating Voronoi regions belong to v a and Vb in a weighted Voronoi 
diagram defined over v a and Vb, with the respective weights d(v a , s),d(vb, s) at v a and Vb- 

Let w(vi), . . . , w(vj),w(vj + i), . . . , w(vk), u>(ffc+i), . . . , w(vi) be successive wavefront segments. 
When an I-cuive(w(vj),w(vj + i)) intersects another l-cuive(w(vk), w(vk+i)) at a point p, due to 
non-crossing property of shortest paths (Lemma 12. ip . there is no need to progress wavefront seg- 
ments w(vj+i), . . . ,w(vk) any further from p. To avoid overlaps of segments within the wave- 
front, we capture I-curve intersection as an event point and update the wavefront by removing 
w(vj+i), . . . ,w(vk) from the wavefront. Therefore, I-curves are also useful in identifying the wave- 
front changes. 

Given that there could be number of boundary cycles, the associations of wavefront segments 
with the boundary edges on every cycle is of interest in determining the interactions. We show the 
following useful property: 

Lemma 2.2 There exists an association A such that the sequence of boundary edges on a boundary 
cycle that are associated with a segment is a contiguous sequence. This is known as the contiguity 
property for wavefront segments. 

Proof: Consider a sequence of bounding edges ej, . . . , ej, . . . , ef., . . . , ej of a boundary cycle. Let 
two sections of boundary ej, . . . , ej and e^, . . . , e\ are associated with wavefront segment w(v a ) and 
another section of boundary ej, . . . , is associated with a segment w(vb). See Fig. [H Also, let the 
shortest path from Vb to s does not pass through v a , and vice versa. 

^From the definition of associations, for every bounding edge e' G {ej, . . . , ej}, there exists 
a point p' on e' which has a shortest path to s via v a . Similarly, for every bounding edge 
e" £ {efc, . . . ,ej}, there exists a point p" on e" which has a shortest path to s via v a ; and for 
every bounding edge e z £ {ej, . . . , e^}, there exists a point p z on e z which has a shortest path 
to s via Vb- But the line segment joining p z to i>&, intersects either the line segment joining p' 
to v a or the line segment joining p" to v a . This is because the point p z occurs on the section of 
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Figure 7: Associations of wavefront segments with convex chains and enter/exit boundaries 

boundary between points p 1 and p" . Suppose the line segment p z Vb intersects the line segment p"v a 
at a point p (argument for the other case is symmetric). Because of the non-crossing property of 
shortest paths (Lemma l2.ip . the intersection point p has at least two shortest paths to s: one via 
v a , and, another via v^,. In turn, there exists two shortest paths from p" (resp. p z ) to s: one via v a , 
and the other via Vb- Hence w(vb) can be associated with e", and, w(v a ) can be associated with e z . □ 

Fig. [7] shows the associations between a section of wavefront W and a section of boundary B 
with I-curves between adjacent segments of W . 

Any corridor enter /exit boundary e is considered as traversed from the first time it got struck 
by a section of wavefront, and from there on e does not participate in associations. In other words, 
each corridor enter /exit boundary is struck by the wavefront at most only once. When a sec- 
tion of wavefront W strikes an untraversed junction J or corridor C, the updates to associations 
of W are termed as wavefront split; whereas the wavefront merger procedure updates the associ- 
ations of sections of wavefront when an edge of already traversed junction J or corridor C is struck. 

The following example describes both the wavefront split and merger. Consider a junction 
J = (ei, e2, e%) in which no edge is traversed yet. Suppose a section of wavefront SW\ struck edge 
e\ when the wavefront is W(d). Then we need to delete e\ from all of its associations, and the 
associations of SW\ needs to be updated to reflect the associations of wavefront segments in SW\ 
to bounding edges corresponding to ei and 63. As shown in Fig. E|a), the section of wavefront 
SW\ is split into section of wavefront SW[ that is associated with e2, and the section of wavefront 
SW{' that is associated with e 3 , for SW{ U SW{' = SW X . 

For d! > d, when the wavefront is W(d'), suppose ei is struck by SW[ as shown in Fig. [He). 
Now based on SW[ I-curves, the section of boundary associated with SW[ comprise C[ U 64 U C'{. 
For d" > a" , when the wavefront is W(d"), suppose e^ is struck by a section of wavefront SW2 (see 
Fig. EIc)). This causes a boundary split. Both the edges e\ and e2 are not considered as existing 
any more as they were respectively struck by SW\ and SW{ in the past. Then based on SW2 
I-curves, we know that C[ and are reachable from SW2. Since there is a section of boundary 
that is common to both SW{ and SW2, we merge these two sections of wavefront. We merge two 
sections of wavefront whenever there is a boundary split. The wavefront merge procedure computes 
the association of C[ U 64 U C" with bunches in SW[ U SW2, based on the proximity of each of 
C{,&4,C{' with SW[ and SW2. Note that each section of wavefront shown in these pictures may 
comprise bunches which are not necessarily contiguous along the wavefront. This is true even for 
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(a) SWi struck ei 



(b) SW[ struck e2 and SW2 struck a 



Figure 8: Wavefront split and merger 

the resultant section of wavefront that is associated with C[ U 64 U C" . 

Property 2.3 A waveform merger is required whenever there is a boundary split. 

We associate either a section of boundary with a segment, or a section of wavefront with a 
g £ dB, such that neither a g nor a w(v) participates in more than one association. This limits the 
number of interactions between the wavefront and dB to 0(n + m). 

Bunches 

The shortest distance computations and updates of associations are proportional to the number 
of associations involved. Given that n could be much larger than m, it would be interesting 
to explore whether the number of interactions/associations is a function of number of obstacles. 
Also, 0(n) wavefront segments could interact between themselves. To reduce the number of such 
possible events, either we need to reduce the number of wavefront segments initiated or show that 
the number of segment intersections are limited to o(n 2 ). The latter is possible due to the following 
structure: 

Definition 2.2 Let Vj,Vj+i, . . . ,Vk be a maximal list of successive vertices along a corridor convex 
chain CC such that for every two neighboring vertices p,q in this ordered list d(s,q) = d(s,p) + 
d{p, q). Then, a bunch B{vj,Vk) is defined as the sequence of wavefront segments w(vj),w(vj+i), . . . , w{vk)- 

The I-curves among segments within a bunch are straight-lines. Two I-curves are said to diverge 
if they are not intersecting, and will not intersect in future as the wavefront expands. 

Lemma 2.3 The intra-bunch I-curves among the segments within a bunch diverge. 

Proof: Consider a bunch B with wavefront segments w(vj), . . . ,w(vi). Suppose Vj lies on corridor 
convex chain C. See Fig. [9l Let p be a point interior to the hull formed by the vertices of C, 
and Vj, . . . ,vi occur in clockwise (resp. counter-clockwise) direction w.r.t. p. Since any two non- 
parallel lines can intersect at one point at most, and any edge v^Vk+i for j < k < I — 1 extended 
in coutner-clockwise (resp. clockwise) direction w.r.t. p intersects all lines induced by edges in the 
set {vjVj+i, . . . ,Vk-iVk}, the I-curves will not intersect in future on wavefront expansion. □ 
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. I-cuve(w(^ ), w(^ +1 )) 



I-curves 



Figure 9: Diverging I-curves 



Let CC, CC be two corridor convex chains. And, let the vertices of CC be {v\, V2, • • • , t>«, i>i+i, 
. . . , Vj = v z ,Vj + i, . . . , Vk, ■ ■ ■ , v n >}. Let v' zout be a vertex of CC . Suppose w(v' zout ) strikes CC at 
a vertex v z along the tangent from v' zout to v z . The tangent v' zout v z defines a sequence S of CC 
vertices, Vj+i, ■ ■ ■ ,v n >, none of which are visible from v' zout . We initiate a bunch B(y z ,v n t), whose 
wavefront segment centers are Vj,Vj+i, . . . ,v n ', and insert B(v z ,v n r) into the wavefront. 

The initial wavefront consists of a (degenerate) bunch with one segment, w(s). Bunches are 
initiated whenever a wavefront segment strikes a corridor convex chain tangentially. As shown in 
Fig. [TUl when a segment from a bunch initiated from a convex chain CC" strikes another convex 
chain CC with v z v' zin as a common tangent, we initiate a bunch B'(v' zin ,v') from v' zin . Similarly, 
when a segment w{v zout ) of a bunch B(v' zin ,v' p ) strikes with v' zout v z as a common tangent between 
CC and CC', a new bunch B(v z ,v n i) is initiated. This process continues until we reach the corridor 
in which t resides. 

We consider a closed corridor as a special case of open corridor. See Fig. [H Let vi,v% be the 
apex points of two funnels. Suppose 1)2 is not yet struck. After the wavefront strikes v%, (two) 
bunches are initiated from v% whenever the wavefront expands from v\ a distance that equals the 
shortest distance from v\ to v%. This distance is precomputed for each closed corridor. The bunches 
initiated from V2 correspond to two corridor convex chains that originate from V2 and define the 
funnel with apex v<i- Let {V2, u\, U2, • ■ • , Uk} and {V2, u^, u' 2 , . . . , «J} be the sequence of vertices 
along the two convex chains of funnel with apex V2- Then the two bunches that are initiated com- 
prise the wavefront segments itf(i>2), w(u%),w(u2), ■ ■ ■ w(uk) and w(ui), w^u^) ■ ■ ■ w{u[) i.e., with the 
wavefront segment originating at V2 is included in only one bunch. 

Hence at any stage of the algorithm, the wavefront is formed by a set of bunches. 

Property 2.4 The wavefront W(d), for any d, is composed of a set of bunches. 

At any stage of the wavefront progression, a bunch B(v z ,v n i) consists of a sequence of segments 
that are already initiated, followed by the rest of the uninitiated wavefront segments (if there are 
any) that possibly belong to B(v z ,v n >) in future. However, for any vertex v in {v z , . . . , v n >}, w(v) is 
part of B(v z ,v n >) if and only if the wavefront strikes v when the wavefront is at Euclidean distance 
d(s,v z ) + d(v z ,v) from s. Each bunch is maintained in a balanced tree data structure. 

As Fig. [TITa) shows, when v z is struck by the wavefront such that the line joining v z and v' zout is 
a tangent to UH, we initiate a bunch B(v z ,v n i) with w(v z ). With the wavefront expansion, when 
w(v z ) strikes v z+ \, we initiate w{v z+ i) in bunch B{v z ,v n >). See Fig. [TlT b). Similarly, w{v z+ 2) 
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Figure 10: Bunch initiation 



is initiated in bunch B(v z ,v n i) with further wavefront expansion. See Fig. [TOc). Although the 
bunch B(v z ,v n i) is inserted into the wavefront when w(v z ) is initiated, the rest of the wavefront 
segments are in B(v z ,v n i) when the wavefront strikes their corresponding centers. The satellite 
data associated with the nodes of the data structure corresponding to every bunch B facilitates in 
determining which wavefront segments within B are initiated. 

Bunch I-curves and Associations 

Since segments within a bunch does not interact between themselves, only the interactions be- 
tween bunches, and the interaction of bunches with dB are of interest. Given that both the number 
of bunches and \dB\ are a function of number of corridors, which is again a function of number 
of obstacles, the total number of events is a function of m. Similar to associations of segments 
with the boundary sections, we are interested in associations between bunches and the boundary 
sections. 

Consider a corridor convex chain or enter/exit boundary g in a boundary cycle. Let B g be the 
set of bunches from each of which g is reachable. Every bunch S G B g is associated with g (or, 
g is associated with bunch S) if and only if a point on g has shortest Euclidean distance to s via 
the center of a segment w € S. The association is defined by the relation: A C Q x B, where Q is 
the set of corridor convex chains or enter/exit boundaries in V, and B is the set of bunches formed 
during the course of algorithm. Again, note that it is not necessary for any two bunches in B g to 
be contiguous in the wavefront. The example in Fig. [T2lfa) shows an association of bunches with 
dB, and how the associations are updated in Fig. [TZTb). 

Lemma 2.4 There exists an association A such that the sequence of boundary edges on a boundary 
cycle that are associated with a bunch is a contiguous sequence. This is known as contiguity property 
for bunches. 
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zout 




bunch 




bunch 



(a) Both the bunch and w(v z ) are initiated 

at v z 



(b) w(v z+ i) is initiated in the bunch 




bunch 



(c) w(v z +2 is initiated in the bunch 
Figure 11: Initiation of wavefront segments within a bunch 




(a) Associations of Bi,w(s), and B2 
are respectively {61}, {&i,&2}, {62} 



tangent 
bunch I-curve 
reachable edge 



(b) Associations of Bi, w(s), and B2 

are respectively {C 2 ,b 3 }, {b 3 ,C 2 \C 3 }, {64, C£} 



Figure 12: Associations of bunches with dB 
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Figure 13: Open corridor with bunches 
The proof for the above Lemma is similar to Lemma [2.21 

To help in associating bunches with dB, we define Voronoi regions corresponding to bunches. 
Let w(vi),w(vi + i), . . . , w(vk) be the segments of a bunch B(v{, Vf.). Let w(v q ),w(v q+ i), . . . , w(v r ) be 
the segments of a bunch B{y q , v r ). Let the bunches B(vi,Vk), B(v q , v r ) be adjacent along the wave- 
front such that the wavefront segments w{yk) and w(v q ) are adjacent. Then the I-curve(-Ufc, v q ) is the 
inter-inter-bunch I-curve between B{vi,Vk),B{v q ,v r ), denoted as l-curve(B(vi,Vk), B(v q ,v r )). 
Although the intra-bunch I-curves are diverging rays (Lemma I2.3P , inter-bunch I-curves could be 
higher-order curves. 

An example in Fig. [13] shows bunches, I-curves among wavefront segments in a bunch, inter- 
bunch I-curves. The following two definitions help in maintaining the associations: 

Definition 2.3 Let B be the set of bunches associated either with a corridor convex chain or 
corridor enter/exit boundary C, i.e., V6 E B, (C,b) € A. When \B\ > 1, a waveform-section for 
C, denoted by WS(C), is the sequence of bunches in B. Note that the bunches in WS(C) need not 
be contiguous in the wavefront. 

Definition 2.4 Let a bunch b be associated to a set S, where each element of S is a corridor 
convex chain or enter/exit boundary. When \S\ > 1, a boundary-section for bunch b, denoted 
with BS(b), is the contiguous sequence of corridor convex chains or enter /exit boundaries in S. 

Consider a corridor convex chain or enter /exit boundary C for which a waveform-section, 
WS{C), is defined. The above partitioning ensures that if a bunch b, which has not yet struck 
C before, strikes a point p € C before striking any other bounding edge in dB, then the bunch 
b is guaranteed to be in WS(C). Furthermore, consider a bunch b for which a boundary-section, 
BS(b), is defined. On wavefront expansion, if the bunch b strikes dB, it would do so only by strik- 
ing a point p € BS(b). In general, a section of wavefront represents either a sequence of bunches 
(not necessarily contiguous) in a waveform-section, or the bunch associated to a boundary-section. 
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Figure 14: Associations, boundary-sections and waveform-sections of a section of wavefront 



The -Redefined below is used in initiating/updating boundary-/waveform-sections and it helps 
in the analysis. 

Definition 2.5 Suppose a bunch a r is associated with a contiguous sequence of corridor convex 
chains/exit boundaries, say Bs, i.e. (e,a r ) 6 A,Ve 6 Bs- If there exists at least one corridor 
convex chain or enter /exit boundary C € S such that C is solely associated with bunch a r then 
RV(a r ) is defined as Bs; otherwise, RV(a r ) is <fi. 

For a bunch b, suppose |i£V(6)| = 0. Hence the bunch b is associated solely to at least one 
corridor convex chain or enter/exit boundary, say C. Then the waveform-section of C, WS(C), is 
obtained as the (not necessarily contiguous) sequence of bunches associated with C. For a bunch b, 
we define boundary-section for b, BS(b), only if |flV(6)| > 1. When defined, BS(b) is the sequence 
of corridor convex chains/exit boundaries in RV(b). The example in Fig. [T4l shows RVs together 
with associations, boundary-sections, and waveform-sections. 

Fig. Q3] shows RVs, associations, boundary- and waveform-sections for a section of boundary 
and a section of wavefront. 

As shown in Fig. [3l when the wavefront is initiated the boundary-section BS of w(s) is 6162^3- 
We compute the shortest distance d' between BS and w(s). When w(s) strikes b±, a waveform- 
section WS± comprising the boundary C[, 64, C[ is associated with w(s) is initiated, and b\ is deleted 
from BS. At the next event point i.e., when w(s) strikes 63, a new waveform-section WS2 com- 
prising the boundary C 3 ' , be, C 3 is associated with w(s) is initiated, and 63 is deleted from BS. We 
continue initiating/updating/deleting waveform- and boundary-sections until the wavefront strikes 
t. 

Boundary-sections are primarily useful when the boundary splits. The data structure that 
stores the boundary sections' facilitates in finding the convex chain vertices that are visible from 
the vertices of wavefront segments. This in turn helps in initiating bunches. 
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3 Algorithm Outline 



Initially, the wavefront is a segment w(s) with radius e, and the boundary-section (also, the initial 
boundary cycles) is the triangle in which s resides. The event points are based on the interaction of 
wavefront with the boundary cycles. The events occur as the wavefront progresses and are main- 
tained in a min-heap, with the corresponding shortest distance at which an event occurs as the key. 
Based on the type of event, an event handler procedure is invoked to handle the event. The event 
points are categorized into the following types: 

Type-I Occurs when a segment w(v) in a waveform-section WS strikes the asso- 
ciated corridor convex chain or enter/exit bounding edge C. This event 
is determined by computing the shortest distance between WS and C. 
The event handler accomplishes the following: 

1. Updating relevant waveform-sections and their associations. 

2. Updating relevant boundary-sections and their associations. 

3. Handling relevant boundary splits if there are any. 

4. Using the new associations, computing the shortest distances and 
pushing corresponding event points to min-heap. 

5. If the event is occurred due to wavefront strike of an enter /exit 
bounding edge of a corridor C , then Type-Ill events are pushed 
to the min-heap: these event points correspond to the wavefront 
progression at which new bunches from the convex chains of C are 
initiated. 

Type-II Occurs when a segment w(v) in a bunch B strikes the associated 
boundary-section BS. This event is determined by computing the short- 
est distance between BS and w(v). The handling procedure is same as 
the Type-I handler. 

Type-Ill Occurs when a segment w(v) in a section of wavefront SW strikes a cor- 
ridor convex chain C at a point of tangency p such that the line segment 
vp is a tangent from v to C. Either a Type-I or a Type-II event could 
cause Type-Ill event. The handler procedure does the following: 

1. A new bunch B is initiated from p. 

2. The associations of bunch B are determined and the relevant 
waveform-sections are updated. 

3. Also, shortest distance corresponding to new associations are 
pushed to min-heap. 
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a i a 2 a 3 a 4 

(a) A bunch (b) Corresponding BHT 

Figure 15: A bunch and its BHT 

Type- IV Occurs when upper convex hull boundary approximations UHi,UHj of 
any two bunches bi,bj within a waveform-section WS intersect for the 
first time. This is determined by computing the shortest distance be- 
tween pairs of bunches within WS. As there could be O(m) bunches 
within a WS, this could lead to 0(m 2 ) computations. However, the 
number of events are reduced by utilizing the special structure in saving 
the waveform-sections (detailed in Subsection I4.2.1| ). The event handler 
primarily does the following: 

1. Suppose SW is the section of wavefront between b, L and bj such 
that t is not located in the region enclosed by SW U {bi} U {bj}. 
Since no wavefront segment in SW can cause a shortest path to t 
without crossing some wavefront segment not in SW, utilizing the 
non-crossing property of shortest paths (Lemma 12. ip . bunches in 
SW are deleted from the wavefront. 

2. Updating relevant waveform- and boundary-sections. 

3. Pushing new shortest distances with the updated waveform- and 
boundary-sections. 

The detailed descriptions of the determination and handling procedures are given in Section [9j 
We continue processing the events scheduled from the min-heap till the sink t is struck. When this 
happens, we compute the shortest path and distance from s to t. 

4 Data Structures 

This Section describes all the required data structures. 
4.1 Bunch Hull Tree (BHT) 

Each bunch is stored in a balanced tree structure, termed as a bunch hull tree (BHT). Suppose 
that there is no bunch originated from a corridor convex chain CC. Let a vertex Vj of CC be struck 
by a wavefront segment w{v' zout ) such that the line v' zout Vj is a tangent to CC. Also, let Vj G W(d). 
As shown in Fig. [TU1 a new bunch B(vj,v n >) is initiated i.e., a bunch hull tree BHT(vj,v n >) corre- 
sponding to B(vj,v n r) is constructed. A wavefront segment w(v) in a bunch is termed as initiated 
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if d(s,v) is determined. Since the wavefront struck Vj and hence d(vj,s) is known, w(vj) is initi- 
ated. At a given W(d), the wavefront segments of a bunch that are initiated are termed as valid 
segments and the other wavefront segments of bunch are termed as invalid segments (see Fig. PTT|) , 
The invalid segments of B(v z , v n i) may possibly become valid in future as the wavefront progresses. 

The leaves of BHT(vj,v n i) consists of w(vj), . . . , w(v n i), wherein w{vj) is the only valid segment 
and no wavefront segment initiated from the remaining vertices. An example is shown in Fig. [15] 
Each tree node p represents an upper hull of the wavefront segments at the leaves of the subtree 
at p. Let q, r be the child nodes of an internal node p of BHT(vj,v n >). Also, let UH q , UH r be the 
upper hulls at q and r respectively. The node p stores the common tangent segment br p between 
UH q and UH r . This common tangent is known as a bridge between UH q and UH r . The upper 
hull at a node p of the tree can be determined from UH q , br p , and UH r . 

When initiating a BHT, only one wavefront segment within that BHT is initiated. Consider 
a leaf node l Vk of BHT{vj,v n i) that corresponds to w{vk), where vt is a vertex of CC and belongs 
to {vj, . . . ,v n i}. ^,From the definition of a bunch, w(vk) is added to BHT(vj,v n i) if and only if 
d(s, Vk) = d(vk,Vk-i) + d(vk-i, + • • • + d(vj + i,Vj) + d{y j, s). Suppose that this is the case. To 
facilitate in automatic implicit insertion of w{vk) to BHT, the data member wpupdate associated 
with l v , stores the negated Euclidean distance along CC between Vj and Vk- In other words, for 
w(vk) to be inserted to BHT(vj,v n i), the wavefront segments at w(vj + i), . . . ,w(vk-i) needed to 
be initiated in that order and w{vk-\) requires to strike Vf.. The root node r stores d(s,Vj) in a 
data member, shortestdist. 

Property 4.1 Let W(d) be the wavefront that caused a bunch B(vj,v n i) to be initiated. Let Vk € 
{vj, . . . ,v n /}. Suppose w(vk) is struck by w(vk~i) £ B(vj,v n /) when the wavefront is W(d'). Also, 
let r be the root node and l Vk corresponds to w{vk) in BHT(vj,v n >). The wavefront segment w(vk) 
in B(vj,v n i) is valid if and only if (d' — r. shortestdist + l(vk) .wpupdate) > 0. 

The root node of BHT(vj,v n i) refers vertex v' zout in tangentstart. 

Lemma 4.1 The variables stored at roots and leaves of all the BHTs are sufficient to compute the 
shortest distance and shortest path from any valid segment w(vk) to s. 

Proof: Let the leaf node l Vk of BHT(vj,v n i) corresponds to Vk- Let W(d) be the wavefront that 
caused the initiation of bunch B(vj,v n i). Suppose w(vk) is struck by w(vk-i) when the wavefront 
is W(d'). Then d' is the d(s,Vk). The shortest path from a point p on w(vk) to s consists of 
line segment pv^, edges along to v^-i, f r -i to v r -2, • • • , Vj + \ to v z , and the tangent from v z 
to root. tangentstart, including the shortest paths computed from v zou f in similar fashion, until 
reaching s. □ 

Each internal node t v stores the maximum wpupdate of its children. Since each leaf node l Vk 
stores the negated Euclidean distance along CC between Vj and v^, a negative wpupdate at t v 
indicates that all wavefront segments stored in the leaves of the subtree rooted at t v are invalid. 

The important data members of a BHT(vj,v n i) are recapitulated w.r.t. a leaf node l Vk , an 
internal node t v , and the root r: 

l Vk .wpupdate Negated Euclidean distance along CC between Vj and 



t v . wpupdate 
r. tangentstart 



maximum wpupdate of its children 
refers to v' zout , given that w{v' zout ) is the wavefront seg- 
ment that struck Vj to cause B(vj,v n i) 
stores d(s, Vj) 



r. shortestdist 
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Initialization 



Let a vertex v z of CC is struck by a wavefront segment w{v' zout ) such that the line v' zout v z is a 
tangent to CC. The following cases needs to be considered: 

(1) Neither w(v z ) nor any w(vj~) from CC is initiated. 

(2) w(v z ) is not initiated but there exists a w(vk) for some € {v z +i, . . . ,v n /}. 

(3) w(v z ) is a valid segment in some BHT. 

(4) w(v z ) is an invalid segment in some BHT. 

For Case (1) (see Fig. [ToTa) ) . we create a new BHT(v z ,v n i). For the Case (2) (see Fig. [ToTb)). 
we neither need to create nor modify any BHTs. Let SW' be the section of wavefront that struck 
Vk, which in turn caused a bunch. Since the sink vertex t is in its own corridor, progressing SW 
and/or a bunch initiated from v z cannot reach this corridor without crossing SW'; but, the non- 
crossing nature of shortest paths (Lemma 12. ip removes this case. We do nothing for the Case 
(3) (see Fig. [T6lfc)); the argument is same as Case (2). For the Case (4) (see Fig. [ToTdl). let 
v z belongs to BHT(vi,Vk), whose leaves are Zj, ■ ■ ■ , l z =lj, ■ ■ ■ ? h- We split BHT(vi,Vk) 
into two BHTs, BHT(vi,v z -\) and BHT(v z ,Vk), using the split procedure mentioned below. 
However, since B{vi, vp.) cannot reach corridor in which t resides without crossing v' zout v z , we delete 
BHT(vi,Vk)- In this case, the split procedure essentially needs to create only BHT(v z ,Vk) from 
BHT( Vi ,v k ). 

Splitting 

Consider a bunch B(vj,v w ) with wavefront segments w(vj), . . . , w(lk), ■ ■ • , w(l w ) in that or- 
der. The process of forming two bunches B(vj, f/e-i), B(vk, v w ) whose wavefront segments are 
w(vj), . . . , w(lk-i) and w(vk), ■ ■ ■ , w(l w ) respectively, is termed as a bunch split. This is required 
due to either of the following reasons: 

(a) Case (4) mentioned with the bunch initialization. 

(b) All the wavefront segments in bunch B(vj,v w ) are no more associated with the same corridor 
convex chain or enter /exit boundary. 

The former is discussed in Case (4) part of bunch initialization. 

Consider the latter. Suppose we require the procedure to form two BHTs: one with leaves 
lj, . . . , lk-i, and the other with leaves l)., . . . ,l w i.e., we wish to form BHT(vj,Vk~i) and BHT(vk,v w ) 
from BHT(vj,v w ). Let p a be the least common ancestor (LCA) of nodes lj, . . . , l^-i, and let pi, be 
the LCA of nodes 1^, . . . ,l w . Then p a with its subtree is the BHT{lj, lk—i), and pb with its subtree 
is the BHT{lk,l w ). For both of these BHTs, we update the wpupdate data member for each 
internal node occurring along the leftmost and rightmost branches i.e., along the split paths. As 
wpupdate of no leaf node is changed, wpupdate of no other internal node needs to be changed. For 
BHT(lj,lk-i), the tangentstart and shortestdist members are same as the root of BHT(vj,v w ). 
For BHT(lk,l w ), tangentstart of root refers to v z with a special flag enabled to denote that this 
bunch is formed during bunch splits which helps in finding shortest paths passing through vt; and, 
shortestdist stores the shortest Euclidean distance along the boundary of CC from Vj to Vk added 
with d(s, Vj). 

Property 4.2 At any W(d), any vertex vi is a leaf in at most one BHT. 
Lemma 4.2 The bunch initiation takes 0{n) time. 

Proof: There are at most O(n) leaf nodes in a BHT. At every node, we spend 0(1) time in 
initializing data members. Hence, the Case (1) takes 0{n). There is nothing to do in Cases (2) 
and (3). As the Lemma 14.31 shows. Case (4) takes O(lgn) time. □ 
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(b) Case (2) 



SW 




(c) Case (3) (d) Case (4) 

Figure 16: BHT Initialization 
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Lemma 4.3 The split operation takes 0(\gn) time. 

Proof: The number of leaves in a bunch are 0(n). Computing LCAs, p a and pb, take 0(1) time. 
As BHTs are balanced trees, updating wpupdate along the split paths take O(lgn) time. Updates 
required at the root node take 0(1) time. □ 



Lemma 4.4 The total number of bunches at any point of execution of the algorithm are 0(m). 

Proof: In Case (1), we create a new bunch. In Cases (2) and (3), we are not creating any new 
bunches. In Case (4), W6 are deleting T id irnmed.ia.tely after creating T new . Therefore, at any 
instance at most two bunches corresponding to a corridor convex chain are alive, one moving to- 
wards/across one corridor enter /exit boundary, and the second towards/across the other corridor 
enter /exit boundary. Since there are 0(m) corridors with each having at most two convex chains, 
there can be at most 0(m) bunches at any instance during the entire algorithm. A split at a junc- 
tion/corridor causes a bunch to divide into (at most) three bunches; once an enter/exit boundary 
b caused a bunch to split, the same b cannot cause split of any other bunch, causing 0(m) splits in 
total. □ 



4.2 Waveform- and Boundary- Section Trees 

The computation of event points require computing the shortest distance for the wavefront to 
strike the boundary dB either at a point on the corridor convex chain or at a corridor entry/exit 
boundary. The Euclidean distances need only be computed between: 

(i) Waveform-section WS and the associated corridor convex chain or corridor enter /exit bound- 
ary C. 

(ii) Boundary-section BS and the associated bunch B. 

To compute these distances, a naive approach would in case (i) compares the distance of each 
bunch in WS with C, in case (ii) compares the distance of each corridor convex chain or corridor 
enter/exit boundary in BS with B. This would lead to a time complexity which is quadratic in 
number of corridors. In order to improve the time complexity, in case (i), we consider a convex hull 
approximation of WS; in case (ii), we consider a convex hull approximation of BS. The convex-hull 
approximation of WS (resp. BS) is defined by the minimum area convex figure enclosing the WS 
(resp. BS). Using the convex representation, we can efficiently determine when the wavefront next 
strikes the boundary. 

Since the convex-hull approximations are to be maintained as the wavefront changes dynami- 
cally, a hierarchical representation of the convex-hull of boundary- and waveform-sections is con- 
structed and maintained dynamically using a balanced tree structure. The data structures to store 
the waveform-sections and the boundary-sections are described in Subsections 14. 2. II and l4~2?2l The 
shortest distance computations using these trees are described in Section [5j 

4.2.1 Waveform Section Tree ( WST) 

The convex hull approximation of a section of wavefront is maintained as a balanced tree. This tree 
is termed as a waveform- section tree, or WST. Each leaf of a WST refers to a BHT. As BHT is 
a convex approximation of a bunch, it is immediate that each leaf node refers to an upper hull. At 
every internal node p of a WST, is maintained an upper hull of the wavefront segments (bunches) 
at the leaves of the subtree rooted at p. Let W(u) be the upper hull (implicitly) stored at a node u 
of a hull tree. Let u be any internal node of WST. Let P(u) be the path from the root to the node 
u. Each internal node u maintains of fset(u) = d — YlwePfu) °f fset(w). Whenever offset(u) < 
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(a) A waveform-section (b) Corresponding WST 



Figure 17: A waveform-section and its WST 

Pk 




Figure 18: Bridge movement with the wavefront expansion 

then W(u) is invalid. See Fig. [TTJ 

As is standard, |14j . at each internal node we will maintain bridges between the two hulls at the 
child nodes. Let q, r be the child nodes of an internal node p of WST. Let S q = • • • , 

and S r = {bj,bj+i, . . . ,bk} be the bunches at the leaf nodes of the subtrees rooted at q and r 
respectively. Formally, we define the bridge br p at node p as a line segment joining two points, pj 
on a valid segment of bj and p^ on a valid segment of such that both the pj and pk are points of 
tangencies. Let LINE(pj,pk) be the line obtained by extending the line segment PjPk infinitely at 
both of its ends. The bridge br p is defined such that the bunches UiG{i k} an belong to one of 
the closed half-planes defined by LINE(pi,Pk). The upper hull UH P at node p is computed from 
the upper hulls UH q , UH r at nodes q and r respectively. The information stored at internal nodes 
is same as in Overmars et al. |14j . 

Let pj lies on wavefront segment w(vj) of bj, and pk lies on wavefront segment w{vt) of bk- The 
slope of the lines VjPj and VkPk are saved at p to facilitate in constructing bridge br p at a given 
wavefront propagation. Consider the maintenance of bridges at internal nodes with the wavefront 
expansion. 

Property 4.3 For pj G w(vj),pk € w(vk), let br p = (pj,Pk) be a bridge at node p of a WST when 
the wavefront is W(d), andpj,pk are the point of tangencies on w(vj) and w(vk) respectively. Let 
p'j be a point at Euclidean distance (\vjPj \ + d! — d) from Vj along the vector VjPj. Also, let p' k be a 
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point at Euclidean distance (\vkPk\ + d' — d) from Vk along the vector VkPk- If v'jiV'k e then 
the line segment p'-p' k is a bridge br' p ofWST at node p when the wavefront is W(d'). See Fig. \18l 

As the wavefront expands, the bridge br p moves as fast as any point on a bunch b r , for i < r < k, 
expands. Therefore, the traversal of pj (or p^) could happen due to some bunch b r , where r < i or 
r > k i.e., b r would be located in a subtree other than the one rooted at p. We handle the bridge 
maintenance in this case by introducing the notion of dirty bridges in Section [TUJ 

Primitive Operations on WST 

In the Overmars et al. [14J structure each leaf node represents a point. However, our data 
structure contains a bunch at the leaf node. Our data structure operations involve insertion of 
bunches, deletion of bunches, tree splitting, and merging of trees. Since all the points defining a 
bunch are contiguous geometrically and do not overlap with other bunches stored in WST, the 
time complexity to insert /delete a BHT from WST is upper bounded by the time complexity to 
insert /delete a point from Overmars structure. Although the bunch itself is maintained as a hull 
tree, note that we never need to traverse BHT while inserting/deleting/splitting/merging hull trees 
in WST. At each node of WST, we maintain the same information as in Overmars et al. |14j 
structure. 



Splitting a convex hull tree is performed on the balanced tree as specified in Preparata et al. 
[2]. The offset at each node on the split path and in the two trees can be computed along the path. 
The combination of two convex hulls is computed by the bridge construction procedure, taking 
offset of the hulls into account. 

Merging two trees is performed by finding the height in the tree at which the two trees are to 
be merged. The offset of the two hulls to be merged at that node can be computed and bridge 
construction follows the procedure in Preparata et al. [2]. 

Lemma 4.5 The upper hull tree of a set of bunches can be maintained dynamically at the worst- 
case cost of 0((lgm)(lgre)) per insertion/ deletion/merge/ split. And, the data structure uses 0(n) 
space. 

Proof: /.From Lemma 14.41 the number of bunches at the leaves of the hull tree at any stage of the 
algorithm are 0(m). The analysis is same as given in Preparata et al. [2], except that we need to 
spend O(lgm) to locate the leaf node I at which we are interested in inserting/deleting a bunch. 
And for constructing bridges at any node along the path from the from I to root, we spend O(lgra) 
time. □ 



4.2.2 Boundary Section Tree (BST) 

The convex hull approximation of a section of boundary is maintained as a balanced tree. This 
tree is termed as a boundary-section tree or BST. Each leaf of a BST refers to an upper hull of 
a corridor convex chain or enter/exit boundary. Every node p of a BST is an upper hull of the 
corridor convex chains or enter /exit boundaries stored at the leaves of the subtree rooted at p. The 
bridges at internal nodes are defined same as in the case of WST. See Fig. [T9j We remember the 
endpoints of bridges as in Overmars et al. [14] . 

In Overmars et al. [H] structure, each leaf node contains a point. However, our data structure 
contains an upper hull at each leaf node (node that enter/exit boundary is a degenerate upper 
hull). Our data structure operations involve insertion and deletion of upper hulls. Since all the 
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(a) A boundary-section (b) Corresponding BST 



Figure 19: A boundary-section and its BST 

points denning any upper hull are contiguous geometrically and do not overlap with other upper 
hulls stored at leaf nodes of a BST, the time complexity to insert /delete an upper hull from BST 
is upper bounded by the time complexity to insert /delete a point from Overmars et al. structure. 
At each node of BST, we maintain the same information as in Overmars et al. structure. 

Lemma 4.6 The upper hull tree of a contiguous set of corridor convex chains or enter/exit bound- 
ing edges can be dynamically maintained at the worst-case cost o/0((lgm)(lgn)) per insertion/deletion/merge/split. 
And, the data structure uses 0(n) space. 

Proof: Since the number of upper corridor convex chains and enter /exit boundaries are 0(m), the 
proof is same as in Lemma 14.51 □ 



4.3 Boundary Cycle List (BCL) 

To help in boundary splits, all the boundary cycles are saved in a data structure. We use a doubly- 
linked circular list, termed as Boundary Cycle List Lbc to denote each boundary cycle BC. Every 
node in Lbc refers to either a BST or WST. Let BS be a contiguous boundary in BC. Let 
ST\ is either a BST or W ST, and ST2 is another BST or WST. Let the sections of contiguous 
boundaries in ST% and ST2 respectively be BS± and BS2- Also, let BSi U BS2 = BS. For every 
two such sections of boundary BS\ and BS2, there exists two nodes in Lbc that are adjacent and 
referring to ST\ and ST2. 

Each event point corresponding to the strike of wavefront with the boundary, pushes to event 
heap a reference to an entry p in list Lbc, s ° that p refers to a BST or WST involved in the 
occurrence of this event. New BSTs or WSTs can be inserted/deleted to this list in O(l) time as 
we refer to the node at which change is happening. Also when boundary splits, to partition one 
BCL into two takes O(l) time only as the computations are local. ^From Lemma 14.41 the space 
requirement is 0(m). 

4.4 Corridors and Junctions 

The partitioned polygonal domain V' into corridors and junctions is modeled as a graph G(V,E). 
The set V of vertices in G comprise the set of junctions in V' , and the set E of edges in G comprise 
the set of corridors in V' . Let e € E be an edge in G corresponding to a corridor C in V' , and let 
v € V be a vertex in G corresponding to a junction J in V' . Then v lies on e if and only if C and 
J share an edge in V' . In other words, two vertices v\,V2 £ V are adjacent in G whenever there is 
a corridor that is adjacent to junctions corresponding to these vertices. Note that the maximum 
degree of G is three. 
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Each edge in G refers to a corridor, which comprise at most two convex chains and two enter/exit 
boundaries. Each of these convex chains is stored in an array data structure, which is envisaged as 
a tree. 

5 Shortest Distance Computations 

5.1 BST and its Association 

The shortest Euclidean distance between a boundary-section BS and its associated bunch B is the 
minimum amount of wavefront expansion required for some segment within the bunch B to strike a 
point in the boundary-section BS. The computation starts at the root node of BST corresponding 
to boundary-section BS. 

Let {li, . . . , lj, . . . , Ik} be the upper hull vertices at node l p , and let {rv, . . . , ry, . . . , r^} be the 
upper hull vertices at node r p so that l p and r p are the left and right children of node p. Also, let 
the line segment l^ry is the bridge at node p. Then we call the hulls formed by the ordered list of 
vertices {li, . . . , lj}, {ry , . . . , ry} as the outer hulls at nodes l p and r p respectively. Similarly, we 
call the hulls formed by the ordered list of vertices {lj, . . . , l^}, {r^, . . . ,rj/} as the inner hulls at 
nodes l p and r p respectively. The definitions of outer and inner hulls are applicable to (implicit) 
hulls stored at any node of BST, other than the root. 

Starting from the root node of BST, we recurse on its nodes as described below. Consider the 
shortest distance computation involved at node p of BST. 

(a) Suppose that bunch B associated with BST does not intersect with the upper hull UH p at 
p. Then the shortest distance from B is to either of the following: 

(1) br p 

(2) outer hull at node l p 

(3) outer hull at node r p 

We compute the shortest distance between the vertices of B and each of these. In Case ([I]), we 
found the required shortest distance. In Case ([!]), we recurse on outer hull at node l p ; in Case (|3|), 
we recurse on outer hull at node r p . The recursion bottoms whenever the shortest distance is to a 
bridge or to a corridor convex chain or to an enter/exit boundary. Suppose recursively traversing 
BST lead the shortest distance to happen at a leaf /, which refers to a corridor convex chain. Since 
we stored each corridor convex chain as a tree, we continue recursing over the tree referred by leaf I 
until we find an edge with which we can compute the shortest distance from B. The case in which 
leaf refers to a corridor enter/exit boundary, is similar. 

(b) Suppose that bunch B associated with BST does intersect with the upper hull U H p at p. 
Then B intersects with either of the following: 

(4) br p 

(5) outer hull at node l p 

(6) outer hull at node r p 

Computing shortest distances from valid segments in B to each of these, we determine which among 
these are traversed. In Case the bridge br p is split to recurse on the inner hulls at l p and r p . 
In Case ([5]), we recurse on the outer hull at l p ; in Case ©, we recurse on the outer hull at r p . 

Now consider the sub-procedure required in (a) and (b): finding the shortest distance d between 
bunch B and a line segment L (or a degenerate point p). To this end, we exploit the unimodal 
property of shortest distance between the wavefront segments in B and L. We do binary search 
over the valid wavefront segments in B to find d, which is similar to computing the shortest distance 
between two convex chains. 
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Lemma 5.1 The shortest distance between the bunch B and a bridge or exit boundary located at 
a node of BST or an edge located at a leaf node of a corridor convex chain (which is a leaf node 
of BST) can be computed in 0((lgm)(lgn)) time amortized over the number of bridges in all the 
BSTs. 

Proof: Since there are 0{m) useful corridors and corridor boundaries together, there are 0{m) 
elements in all the BSTs together at any point of execution of the algorithm. Hence, every BST 
is of O(lgm) depth during the entire algorithm. 

Finding the shortest distance between the bunch B and a bridge/edge takes O(lgn) time by 
doing binary search over the possible O(n) intra-bunch I-curves. 

Suppose that there is no change in both BST and B. Then no case mentioned under (a) and 
(b), takes more than 0((lgn)(lgm)) time, as in the worst case we need to traverse BST along its 
depth doing binary search over the I-curves of B during this traversal. In Case (j3J), we may traverse 
the inner hulls of both the nodes at every stage. But that traversal can be charged to the bridge 
that was intersecting B. Once a bridge is split it won't split again, as each bridge is treated similar 
to an edge being struck by the wavefront. 

Suppose the shortest distance from bunch B is found to a leaf node of BST, which is a corridor 
convex chain, CH. We can compute the shortest distance between BHT and CH in O(lgn) time. 
The exit boundary stored at a leaf node of BST is treated in the same way as a bridge of BST. 

□ 



5.2 WST and its Association 

The shortest Euclidean distance between a waveform-section WS and its associated corridor convex 
chain or enter/exit boundary C is the minimum amount of wavefront expansion required for some 
wavefront segment within the waveform-section WS to strike C. 

Let us consider the case in which WS is associated with a corridor convex chain C. Let the 
tree corresponding to C be HT. Since the intra-bunch I-curves diverge, inter-bunch I-curves in 
the given waveform-section play a major role. The computation starts at the root node of WST 
corresponding to the waveform-section WS. At a node p of WST if corridor convex chain C does 
not intersect the bridge at p, then the shortest distance from C to WS is either to the bridge at p 
or lies to an upper hull at its left/right child subtree root. In the first case, the shortest distance 
is noted as the distance between the bridge and the corridor convex chain C . If the convex chain 
stored at left (resp. right) descendant node is nearer, recurse on the left (resp. right) subtree. In 
the case that the corridor convex chain C intersects the bridge at p of WST, the bridge is split 
to find the shortest distance between the C and the upper hulls located at the subtrees of p. To 
consider the case in which some point q G C lies inside an upper hull ch p stored at node p of WST, 
we denote the distance to a bridge of WST from q as positive (resp. negative) whenever q lies 
outside (resp. inside) ch p . Let br be the closest bridge to C located at a leaf node of WS, and 
this leaf node corresponds to a bunch B with bunch hull tree as BHT. Then, the computation 
recurse on BHT similar to the traversal of WST explained above, considering only bridges over 
the valid segments in bunch B. This is accomplished by noting that for a non-root node p of BHT, 
the negative root.wpupdate + p.wpudate indicates that the rooted subtree with p as root has only 
invalid segments as leaves. 

The other case in which the waveform-section WS is associated with an exit boundary e is same 
as the computation where the corridor convex chain CH is having only one edge. 

Suppose the computed shortest distance for a waveform-section WS to reach its associated 
corridor convex chain or enter /exit boundary C is represented as an event point evt in the event 
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heap. Let two inter-bunch I-curves of waveform-section WS intersect before the evt occurs. This 
require updates to both the WST and the shortest distance, which are attached as satellite data 
associated with evt. However, the intersection of two inter-bunch I-curves results in a Type-IV 
event, which in turn takes care of these updates. 

Lemma 5.2 The shortest distance between the corridor convex chain or enter/exit boundary C 
associated with the waveform-section WS, and a bridge/segment located at a node of WST (or, a 
node of BHT in WST) can be computed with 0((lgm)(lgn)) time complexity amortized over the 
number of bridges in all WSTs. 

Proof: We analyze the case in which WS is associated with a corridor convex chain; the other 
case when WS is associated with an exit boundary is similar to this. According to the definition 
of WST, the waveform is formed with bunches at the leaves. The total number of bunches at any 
point of execution of the algorithm are 0(m) (from Lemma l4.4p . Hence, any WST is of O(lgm) 
depth during the entire algorithm. For a bridge br at node p of WST, we need to compute the 
shortest distance between C and br, between C and the two upper hulls implicitly stored at p's 
immediate child nodes to find the closest bridge. Since WST is of O(lgm) depth and the tree 
corresponding to C is of O(lgn) depth, every shortest distance computation takes 0((lgm)(lgn)). 
This complexity includes the computation required in determining whether the closest point pt € C 
lies inside/outside the upper hull stored at node p. Once a bridge br in WST is determined as 
closest to C , br will never be considered again. 

We next consider the complexity of computing shortest distances' within bunch hull trees. Let 
BHT is a bunch hull tree located at a leaf node pi of W ST. Also, let the bridge at node pi is 
determined as closest to a bridge br' in HT among all the bridges that are not struck in HT. Then 
we need to traverse the bridges over the valid segments of that BHT; similar to the traversal of 
WST bridges, this is of O(lgn) amortized complexity. □ 



6 Ilntersect Procedure 

Let SB = {Si, S2, ■ ■ ■ , Sk} be the sequence of contiguous sections of boundary such that Sj rep- 
resents either a boundary-section associated to a bunch, or a corridor convex chain (or enter/exit 
boundary) associated with a waveform-section. The Ilntersect procedure finds the corridor convex 
chain or enter /exit boundary in SB with which a given I-curve(u;(ci), w(c2)) first intersects. We 
find the intersection of l-cuxve(w(ci), w(p2)) with each of S\,S2, ■■■ ,Sp, in that order (the cost is 
amortized). If the I-curve(i«(ci), wfa)) does not intersect with SB, the procedure returns the 
same. Since an I-curve is not a straight-line, it may intersect a SB multiple times. However, we 
are interested in the first point of intersection along the given I-curve(w(ci), w{c2)) with SB. 

First, consider finding the intersection of I-curve(u;(ci), w{c2)) with Si £ SB. Let HT be the 
hull tree of Sj. We denote the upper hull at a node p of HT with H p , outer upper hull at p with 
OH p , and the inner inner hull at p with IH p . The computation starts at the root node of HT. 
Consider the computation at a node p of HT. Let q, r be the left and right children of p. The bridge 
br p at p be piP2- The following are the possible places at which I-curve(w;(ci), w{c2)) intersects: 

(1) OH r 

(2) OH q 

(3) br p 

(4) None 

^From the definition of I-curve, if I-curve intersects bridge br p then p\ (resp. P2) is closer to 
ci and P2 (resp. p\) is closer to C2- This fact is used in differentiating Case ([3]) from the first 
two Cases. Let the angle c\p\ (resp. C2P1) makes at c\ (resp. C2) is less than the angle made by 
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C\P2 (resp. C2P2) at c\ (resp. C2). The conditions in which Case ([T]) (resp. Case ((2J)) occur is: 
|cipi| < I C2P1 1 , 1 C1P2 1 < IC2P2 1 and the line segment c\p\ (resp. C2P1) is in the left (resp. right) 
half-plane defined by the line C2P2 (resp. C\P2)\ or, \c\p\\ > | C2P1 1 , \cip2\ > \c2P2\ and the line 
segment c\p\ (resp. C2P1) is in the left (resp. right) half-plane defined by the line C2P2 (resp. c\p<z)- 
The Case ([3|) occur if p\ is closer to c\ (resp. C2) than C2 (resp. ci). The Case ^ occur if we reach 
a leaf node and do not find the intersection point. 

For Case ([!]) (resp. Case ([2])), it is guaranteed that the I-curve does not first intersect with any 
of the elements located in the right (resp. left) subtree of r. So we recursively traverse left (resp. 
right) subtree of r. For the Case ([3]), we need to split the bridge br to further find the intersection 
of I-curve with either of IH q ,IH r . Since we do not know which inner hull is intersected by the 
I-curve, we traverse both left and right inner hulls at r. To facilitate traversing inner hulls, we 
determine the visible cones from both the centers c\ and C2 by computing tangents to inner hulls. 

The worst-case occurs when we need to compute the tangents to inner hulls along the depth 
of HT, leading to the time complexity of 0((lgm)(lgm)). Since the I-curve intersects at most 
O(lgm) bridges and no bridge splits more than once, the amortized time complexity is O(lgm). 
Since the number of boundary- and waveform-sections are bonded by 0(m) i.e., since the number 
of elements in SB is 0(m), finding the intersection of the given I-curve with the SB amortized 
over splits, mergers and Type-IV intersections takes 0(lg 2 m) amortized time. 

7 Merging 

Consider a section of boundary SBa associated with a section of wavefront A. Also, consider a 
section of boundary SBb associated with a section of wavefront B. Suppose B struck e, causing 
the reachable edges of B to include new bounding edges from SBa- Then we need to update the 
associations of B by considering the proximity of bunches in A U B with the corridor convex chains 
or enter /exit boundaries in SBa such that: 

1. the sections of boundary associated with both A and B are contiguous, 

2. each element in SBa is associated with either ^4 or B, 

3. at most one element in SBa is associated with both A and B. 

The process of updating the associations of A and B and initiating/updating boundary- and 
waveform-sections for all elements in SBa U SBb is termed as the merger of sections of wave- 
fronts' A and B. 

Consider Fig. [20T a). Here, the section of wavefront A comprising bunches a\,...,a%, and 
SBa comprise corridor convex chains or enter/exit boundaries ci, ... , eg. Specifically, BS{a{) = 
{c±, . . . ,c§}, BS(a,2) = {c5,C6,C7}, and BS{a$) = {cj, cs}- Suppose an event caused part of 
SBa to be reachable from a waveform-section B consisting of bunches b\,..., 65. The asso- 
ciations of SBa U SBb after the wavefront merger is shown in Fig. [201b) with inter-bunch 
I-curves: BSfa) = {ci,c 2 }, BS(b 2 ) = {c 2 ,c 3 ,c 4 }, WS(c4) = {b 2 ,b 3 ,ai}, BS(a 1 ) = {c 4 ,c 5 }, 
BS(ci2) = {cs,cq,cj}, and BSi^a^) = {cj, cs}- Also, the corresponding boundary- and waveform- 
section trees are initiated/updated. As mentioned below, new associations follow all three charac- 
terizations: 

1. the contiguous section of boundary c±, . . . , C4 is associated with B, and the contiguous section 
of boundary C4, . . . , c$ is associated with A, 

2. every element in ci, . . . , eg is associated with either A or B, and 

3. only C4 is associated with both A and B. 

The inter-bunch I-curves, bunches, and the corresponding sections of boundary associated to- 
gether define Voronoi regions. For example, before invoking the merge procedure (Fig. [207a)). the 
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(b) After merge: ci, . . . , ca are associated with A, and C4, . . . , eg are associated with B 



Figure 20: Modification to associations of two sections of wavefronts in Merge procedure 
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Voronoi regions are: I-curve(ao, ai), a\, I-curve(ai, 02), section of boundary along SBa between 
po and pi; I-curve(ai, 02), 02, I-curve(a2, 03), section of boundary along SBa between p\ and p% 
I-curve(a3, 04), 03, I-curve (03, CI4), section of boundary along SBa between P2 and p%. After invok- 
ing the merge procedure (Fig. [20T b)). the Voronoi regions are: I-curve(6o, 61), b±, I-curve(6i, 62), 
section of boundary along SBa between p' and p%; etc.,; I-curve(63, 64), 64,01, I-curve(ai, 02), sec- 
tion of boundary along SBa between p' 3 and p%; etc.,; At most one Voronoi region is bounded by 
one inter-bunch I-curve from A and another from B. In updating the Voronoi regions, some of the 
inter-bunch I-curves from A and B are not used. For example, although I-curve(&4, 65) intersects 
with eg, section of boundary consisting of C5, . . . , eg is closer to A than the bunches in B. Although 
this example shows all the bunches in A (resp. B) are contiguous along the wavefront, it is not a 
requirement. 

Consider any corridor convex chain or enter /exit boundary, say C3 in SBa, that is associated 
with bunch a\. For wavefront segment w(a' 1 ) G a\, let d! be the Euclidean distance between 
w(a'i) and C3. Given the intersection of inter-bunch I-curves of B with SBa, we know that C3 
is closer to bunch 62 among all the bunches in B. For a wavefront segment w{b' 2 ) € 62, let d" be 
the Euclidean distance between w(b' 2 ) and C3. If d" < d' , then the association of C3 is changed to 62. 

Herewith we describe the procedure to associate a section of SBa to bunches in B. For each 
of the corridor convex chain or enter/exit boundary c occurring in SBa, the proximity of c to its 
association in A requires to be compared against the proximity of c with every bunch in B, and 
re-associated when necessary. Before B strike SBa, suppose that there is no Type-IV event due 
to bunches in B. Then the intersection points of inter-bunch I-curves from B are in sorted order 
along SBa- For example, in Fig. [201 I-curve(&3, 64) is intersecting C4; if c±,... ,04 are associated 
with 61,62,63, then the section of boundary c^,...,c% is definitely not associated with 61,62,63. 
This unimodal property accommodates binary search over the inter-bunch I-curves of B in updat- 
ing the associations. Since the inter-bunch I-curves are higher order curves, we use the Ilntersect 
procedure described in Section [6] in determining which element of SBa is intersected by an I-curve 
under consideration. 

The following Lemma says that there exists a contiguous section of boundary that can be 
associated with a section of wavefront. 

Lemma 7.1 There exists an association A such that the sequence of boundary edges on a boundary 
cycle that are associated with a section of wavefront is a contiguous sequence. This is known as the 
contiguity property for sections of wavefront. 

Proof: Immediate from Lemmas 12.11 and 12.41 □ 

This property says that there always exists an implicit I-curve between A and B that separates 
the contiguous section of boundary that is associated with A from the contiguous section of bound- 
ary that is associated with B. Suppose all the elements that occurred before c € SBa along SBa 
are associated with B, and c needs to be associated with both A and B (or, with A only). Then 
no element in SBa that occur after c needs to be associated with B. This is considered in halting 
the binary search procedure. 

Here is the description of ASSOCATOB procedure. Let the bunches in B be b\, . . . , b q . Starting 
from 61, the procedure finds a bunch 6,- in B such that between the points of intersection of inter- 
bunch I-curve ( 6j _ 1, 6j) and the inter-bunch I-curve(6j, 67+1), there exists a contiguous section of 
boundary SB' £ SB A and \SB'\ > 1 (line 4 of ASSOCATOB). If such an I-curve is found, RV(bj) 
is set to SB'; and, RV(b\), . . . ,RV(bj) are set to 4> (lines 5-8). In the next iteration, the procedure 
starts from bj. The iterations terminate when an element in RV{A) needs to be left with its current 
association (line 11). Then it adjusts the RVs of A; and, updates BSTs and WSTs of bunches 
in both A and B (lines 12-13). After the change in associations, the event points corresponding 
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Figure 21: Associating SV A U SV B toAuB 



to the shortest distance between boundary-sections and their associated bunches, and the shortest 
distance between waveform-sections and their associated corridor convex chains/exit boundaries 
are updated in the event min-heap. Also, for every corridor convex chain C whose associations are 
changed within the procedure, similar to the procedure described in Type-I/Type-II events, we find 
the tangent to C and push the corresponding Type-Ill event to min-heap (line 14). Note that the 
Type-Ill event generated herewith could cause Case (4) of BHT initialization (See Subsection l4.ll ). 

Also, suppose C € SB a is the corridor convex chain or enter /exit boundary encountered by 
the binary search procedure in ASSOCATOB, and C is found to be associated with A. Then the 
procedure does not process elements in SBa that occur after C. Let A R (resp. B R ) denote the 
sequence bunches in A (resp. B) ordered in reverse direction. Since the procedure does not know 
from which end of the sequence of bunches in B that it requires to start associating with SBa, we 
need to invoke with arguments: A, B; A,B R ; A R ,B; A R ,B R . Although explicitly we do not state, 
all these eight invocations of ASSOCATOB are considered in MERGE procedure. 

As mentioned earlier, we invoked the merge procedure when B struck an edge e. Suppose B 
was associated with a section of boundary SB' B = SBb U {e}. When it strikes e, the section of 
boundary SBb is reachable from A. See Fig. [21] Then the merge procedure needs to extend the 
three characterizations given in the beginning of this Section to SBa U SBb and A U B, so that 
to associate SBa U SB b with the bunches in A and B. In other words, we also require to invoke 
ASSOCATOB with arguments: B,A; B R , A; B,A R ; B R ,A R . Considering all these possibilities, 
the MERGE procedure invokes ASSOCATOB eight times for merging the two given sections of 
wavefront. 

As a whole, the merger takes two sets, say A and B, of bunches and their associations in terms 
of WSTs and BSTs as input, and outputs a set C C A U B of bunches with initiated/updated 
WSTs and BSTs. 

Theorem 7.1 The ASSOC 'ATOB(A,B) procedure correctly updates the associations of bunches in 
both the sections of wavefronts' A and B. 

Proof: Let SBa and SB B be the sections of boundary associated with bunches in A and B respec- 
tively. We intend to prove three characterizations given in the beginning of this Section. Primarily, 
we prove that for every bunch b in AUB, section of boundary associated with b is contiguous along 
a boundary cycle whenever SBa U SBb is contiguous before invoking the procedure. 

The procedure assumes that both SBa, SBb are contiguous along a boundary cycle. To justify 
this assumption, note that the association of bunches are updated either in ASSOCATOB proce- 
dure or in the procedure which splits a section of wavefront. Since the wavefront split procedure 
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Algorithm 1 Associates corridor convex chains or enter/exit boundaries in RVs of a section of 
wavefront A to RVs of another section of wavefront B; also, updates WSTs and BSTs for the 

bunches in B 

procedure ASSOCATOB(A, B) 

Require: the re-association starts by considering the first bunch in each of A and B. The bunches 
in both A and B are considered sequentially from there on; the procedure terminates when it 
finds the first convex chain or enter /exit boundary among the sequence of convex chains/exit 
boundaries of RV(A) which cannot be associated with a bunch in B 

1: Let 

A : a±, ct2, ■ ■ ■ , a p , a sequence of bunches 
B : b±, 62, ■ ■ ■ , b g , a sequence of bunches 

RV{A) : {Jl=iRV( a i) = c i> C2, . . . , c r , contiguous sequence comprising corridor convex 
chains/exit boundaries in dB (note that this sequence is not computed explicitly) 

2: prevb := 0, prevc := 

3: repeat 

4: by doing binary search over the range [prevb + 1, q] find j, where j = mm{i : i € [prev b+ 1, q] 
and I-curve(6j, satisfies c s := IIntersect(I-curve(6j, RV(A)), c s 7^ c pre vc} and points 
on c s can be associated with bi instead of at where currently c s _i and c s € RV{ak)- 

5: set RV(bpr evb+1 ),. . . ,i?V r (6 j _i),i?y(6 i ) to 0. 

6: if a corridor convex chain or enter/exit boundary c € {c prevc+ i, . . . ,c s } can be entirely 

associated with bj then 
7: add {c prevc+ i, . . . , c s } to 
8: end if 

9: prevb:=j, c prevc :=c s 
10: until prevb > q — 1 

11: remove corridor convex chains or enter/exit boundaries from RV(A) so that RV(A) C\RV(B) = 

4>. 

12: for every bunch bi € B, if RV(b,j) is modified and = 0, and further, bi is associated 

to a corridor convex chain or enter/exit boundary c then update WST(c) s.t. WS(c) is the 
sequence of bunches associated to c. 

13: for every bunch bi G 5, if -Ry(fri) is modified and \RV(pi)\ > 1 update BST(bi) s.t. BS(bi) is 
the same as 

14: for every corridor convex chain C whose associations are changed, determine and push a Type- 
Ill event to min-heap. 



Algorithm 2 Modifies RVs, WSTs and BSTs for the two given sections of the wavefront: 

SW 1 ,SW 2 

procedure MERGE(SWi, SW 2 ) 

l: ASSOCATOB(Wi, SW 2 ) 
2: ASSOCATOB(W 2 , SWi) 
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maintains the contiguity property, the input bunch associations are guaranteed to be contiguous, 
provided that ASSOCATOB procedure maintains the contiguity. It also assumes that SB a U SBb 
is contiguous before invoking the procedure. This is valid since we always invoke the procedure 
with arguments which obey this precondition. 

First, note that there exists a contiguous association of SB SBb to bunches in AL)B, due to 
Lemma [7. 11 The re-association of SBa to B starts from the first element in SB a- Consider the first 
bunch b\ in B. The procedure assigns all the corridor convex chains or enter/exit boundaries prior 
to the intersection of I-curve(&i, 62) with SBa to b\, if these boundary elements are determined 
closer to b\ as compared to segments in A. Let this set be SB'. Since the associations start from 
the first element in SBa, and the elements in SBa U SBb are contiguous, it is guaranteed that 
SB' together with the section of boundary already associated with bunch b\ prior to the invoca- 
tion of this procedure is contiguous. Inductively, this argument can be extended to all the bunch 
re-associations in B. The binary search (line 4) over the inter-bunch I-curves returns the correct 
next I-curve which does not intersect c prevc . Once j in bj is incremented (line 4) it never gets 
decremented, hence associating boundary edges from A to B (line 7) maintains the contiguity for 
each bunch. In every iteration of the repeat loop (line 3), we iterate sequentially along the I-curves 
among bunches in B. Hence SBb is contiguous after the re-associations. All the convex chains or 
enter/exit boundaries removed from SBa (line 11) are contiguous, hence SBb is contiguous after 
the invocation of ASSOCATOB. For the bunches in A whose associations are left as they were, 
based on the preconditions, contiguity is obviously maintained. 

Note that ASSOCATOB procedure does not take into account the case in which two inter-bunch 
I-curves of B intersect before they intersect the section of boundary associated with A. The correct- 
ness for this case relies on the fact that the boundary- and/or waveform-sections of bunches in B are 
updated whenever the Type- IV event corresponding to this inter-bunch I-curve intersection occurs. 

Also, the correctness of Ilntersect procedure is discussed in Section [6l □ 

Corollary 7.1.1 The MERGE(A, B) procedure correctly updates the associations of bunches in 
both the sections of wavefronts' A and B. 

8 Boundary Split 

As described in Section [21 a boundary cycle may split. Suppose the edge e of a junction J (or C) is 
associated with the wavefront from both of its sides. This causes e to appear twice in a boundary 
cycle BC. In Fig. [51 the bounding edges b{ and bj of BC correspond to e. Whenever edge e 
is struck from either of its sides, the boundary cycle splits into two disjoint boundary cycles, say 
BCi and BCj. Let REG \, REG 'j be the regions bounded by BC{ and BCj respectively. Both the 
boundary cycles BCi and BCj are surrounded by the wavefront, whereas at least two edges of J 
(or C) are struck. Considering the non-crossing property of shortest paths (Lemma l2.ip : 

• a shortest path to t could either traverse only some edges / vertices of REGj , or 

• traverse some edges/vertices of REGi, one/two edges of J (or C), or 

• traverse some edges/vertices of REGi, one/two edges of J (or C), and some edges/vertices 
of REGj in that order. 

In other words, the only way in which a shortest path traversing REGi can reach t is through J 
(or C) as shown in Fig. [221 Hence we refer J (resp. C) as a gateway, say g. 

Definition 8.1 For a bounding edge e £ S of a junction/corridor g, suppose e is the last edge that 
struck, which resulted in the split of a boundary cycle BC into BCi an d BCj, then g is termed 
as a gateway. For every edge e' of g, if e' incident to both BCi an d BCj then e' is termed as a 
gateway-edge of g. 
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We denote that the gateway g is attached to REGi (resp. REGj,BCi,BCj), or REGi (resp. 
REGj, BCi, BCj) is attached to g. Taking Euclidean metric with the non-crossing nature of shortest 
paths (Lemma 12. ip into consideration, each gateway is assigned a orientation, so that it suffice to 
expand the wavefront traversing through it over only one region attached to it. Since in our case 
the destination is in REGj, this is the region needed to be traversed by the wavefront expanding 
through g. The gateway g is termed as the outgoing gateway w.r.t. REGi (resp. BCi). 

Lemma 8.1 Every region REG either contains the destination t or is attached with one outgoing 
gateway. 

Proof: We prove this statement using induction. When the algorithm starts, there is only one 
region, which does contain t. Assume that the induction hypothesis holds for the input polygonal 
region having k regions. Now we extend the argument when a boundary split occurs, creating k + 1 
regions. Let a region REG is split into two regions, REGi and REGj. Suppose g' is the gateway 
attached to both REGi and REGj. ^From the induction hypothesis, we know that the region REG 
either contains t or an outgoing gateway g is attached to it. For the former case, where REGj 
(resp. REGi) contains t, we orient g' to be an outgoing gateway w.r.t. REGi (resp. REGj). In 
the case that t is not located in REG and g is attached to REGj (resp. REGi), again we orient g' 
to be an outgoing gateway w.r.t. REGi (resp. REGj). □ 

The wavefront propagation and the initiation of boundary cycles can be represented as a tree. 
The initial shortest path wavefront w(s) is of e radius, and the untraversed region contains t. Due 
to the wavefront progression, suppose REG splits into two, now one of these regions contains t. 
Let us denote the region containing t by REGt and the other region by REG\. Both REGt and 
REG\ are connected with a gateway g%. The orientation of g\ is from REG\ to REGt- Then we 
represent REG\ and REGt as nodes and gateway g\ as an edge in a graph as described below. 

At any point of the algorithm, all the untraversed regions together with the gateways that 
connect respective regions together are represented with an oriented Boundary Cycle Tree, 
BCT{V, A), where the set V comprises of the set of untraversed regions in the polygonal domain and 
the set A comprises of gateways, such that every arc a = (v r , v") € A represents a gateway from the 
boundary cycle represented at v' to the boundary cycle represented at v" . The nodes and edges are 
added to the (logical) boundary cycle tree as the algorithm proceeds i.e., whenever there is boundary 
split, one new node corresponding to the new region and one edge corresponding to the gateway 
are added to BCT. For example, in Fig. [231 consider a path REG t , REG X , REG 12 , REG121 in 
BCT. This corresponds to a boundary splits that occurred over time among the regions along this 



35 



Figure 23: Boundary Cycle Tree and the corresponding region splits 



path, the first split being at the root. Also, a shortest path to t that occurs along this path must 
traverse across a suffix of the regions corresponding to nodes REG121, REG12, REGi, REGt in the 
boundary cycle tree. In other words, consider a section of wavefront W that traverses the region 
at node v of BCT; for W to cause a shortest path SP to t, SP must traverse across the regions 
and gateways at all the ascendant nodes of v. 

A shortest path from s to t may possibly goes through the region corresponding to a node in the 
tree and traverses across all the gateways occurring in the path from that node to the root. At every 
node v of the boundary cycle tree, we merge all the wavefronts that could traverse the gateway 
corresponding to v, say g. In other words, all the wavefronts that traverse the regions/gateways 
associated with the nodes of the subtree rooted at v are merged at g. 

To determine the location of either t or the outgoing gateway w.r.t. a boundary cycle BC, we 
use the algorithm given in |16j . This facilitates in orienting the gateways. 

Lemma 8.2 Given that there exists a path from s to t, a shortest path can be found using some 
( none ) of the gateways at any stage of the wavefront progression. 

Proof: The proof is by induction on the number of regions present when the shortest path from s 
to t is found. Consider the base case in which t resides on or within one boundary cycle by the time 
the shortest path is found. Here, we find a shortest path from s to t without using any gateways. 
Assume that the induction hypothesis holds for the given polygonal region having k regions. Now 
we extend the argument for k + 1 regions. Let a region REG be split into two regions, say REGi 
and REGj. Suppose the region REGj has either t or an outgoing gateway g' w.r.t. REG attached 
to it. Also, suppose g is the outgoing gateway w.r.t. REGi attached to both REGi and REGj. 
See Fig. [M 

^From the induction hypothesis, we know that with k gateways, there is a shortest path from 
s to t where that shortest path uses some (none) of the gateways. If this shortest path does not 
traverse across REG, there is nothing to prove. Otherwise, this shortest path may go in one of 
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Figure 24: A path in the boundary cycle tree 



these ways: traversing across REGj only; traversing across REGi first and then entering REGj; 
traversing across REGj first and later entering REGi before reentering REGj. In the first case, 
we can find a shortest path without using g. For the second case, the orientation of g with the 
non-crossing property of the wavefront assure that such a shortest path is retained. The Euclidean 
metric makes the last case inessential, hence the orientation of g eliminates this case altogether. □ 

Let A and B be the sections of wavefronts associated with the edges bi and bj (the bounding edge 
that is just struck) respectively. Based on the above Lemma, A and B need to be associated with 
BCj only. Let W be the set of bunches associated with boundary cycle BCj that either contains 
t or attached with an outgoing gateway. Let A' denote the listing of bunches in W ordered so as 
to end with the last bunch A associated with the bounding edge bi while traversing the wavefront. 
Similarly, we define B' to be the wavefront bunches ordered in the reverse direction so as to end 
with bj. We merge A' and B' by invoking the MERGE procedure. This determines the association 
of bunches in A and B to the bounding edges of the boundary cycle BCj. 



Wavefront propagation along the gateways 

We next show how to process the wavefront propagation over the boundary cycles. The orien- 
tation of the gateways imposes an ordering of boundary cycles for processing the regions. Consider 
the current wavefront associated with region REGi and the outgoing gateway g attached to it. We 
expand the shortest path wavefront over REGi, until every corridor /junction bounding edge in that 
region is struck. When this happens, we say that REGi is traversed. When BCi is a degenerate 
cycle comprising one vertex, which is not t, this traversal is not required and the algorithm can 
proceed to processing of REGj . During the traversal of REGi , all the bunches that strike gateway- 
edges of g are combined into one section of wavefront, termed as B. Suppose g' is the outgoing 
gateway attached to REGj. Also, let g' i,g' 2 be the sides of edges of g from outside of g. Consider 
any region REG' corresponding to an ancestor of node corresponding to REGi in the boundary 
cycle tree. See Fig. [MJ ^From non-crossing property of sections of wavefront (Lemma 12. ip . note 
that there is no way for a wavefront segment in B to reach a point located in REG', without 
traversing some of: g' ol , g' o2 , edges in BCj. Since wavefront segments in B could potentially expand 
and traverse into REGj and the regions that are ancestors of REGj, the section of wavefront B is 
associated with the bounding edges of BCj,g' ol and g' o2 so that to cover all possible non-crossing 
shortest paths. 

We backtrack in time and restart expanding the shortest path wavefront from the time at which 
BC was split into BCi and BCj, say time tj. During this re-traversal, we are interested in associa- 
tions of B with REGj , g' ol and g' o2 . Let A be the sections of wavefront associated to the bounding 
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edges of BCj, g' ol and g' o2 . We invoke MERGE procedure to merge B with A. Let A' denote 
the listing of bunches in A starting with the first bunch among the bunches in A associated with 
the bounding edge Then we invoke MERGE procedure with A' and B as arguments. The 

combined wavefront will then expand into the region REGj, starting at time tj. 

However, note that B is the result of wavefront propagation that specifies the bunches at a 
future time t', whereas A is the wavefront section when the boundary split occurs (i.e., at time tj). 
We thus create an offset at the root of B equal to time —t'. This negative value ensures that the 
wavefront section B and A are processed at the same time. A segment that has a negative offset 
is not active and will not be considered for the shortest distance computations. The computations 
will take into account the offsets of segments in B. 

Lemma 8.3 The associations of all sections of wavefront involved in an invocation of a merge 
procedure are updated correctly. 

Proof: Suppose the boundary cycle BC is split into BCi and BCj. Let REGi, REGj be the regions 
bounded by BCi, BCj, respectively. Let g',g be outgoing gateways w.r.t. BC,BCi respectively. 
Further, let g' \,g' 2 be the sides of g' from outside g' . The region REGj is guaranteed to either 
contain t or attached with an outgoing gateway. In both these cases, from Lemma 18.21 it suffice to 
associate the sections of wavefront traversing through gateway g solely to the bounding edges of 
BCj, g' ol , g' o2 . Hence, while progressing the wavefront over REGi, we combine all the sections of 
wavefront associated with the gateway-edges of g into a section of wavefront. And, we associate 
this section of wavefront with the bounding edges of boundary cycle BCj and the gateway-edges of 
g' . To reduce the number of event points, at first we traverse REGi and then we determine these 
associations by invoking MERGE procedure. After these associations are determined, we offset the 
wavefront segments that traverse region REGi by the time difference from the time at which BC 
was split. 

Suppose a bunch w was determined to be in a section of wavefront B which struck gateway- 
edges of g. As the wavefront progresses, the bunch w may get split or merged with other sections 
of wavefronts. Since we are, in effect, restarting the wavefront from the time at which BC was split 
into BCi and BCj, the correctness proof needs to take into account the following: 

1. bunch w may either not be alive. 

2. bunch w associated with an untraversed edge in BCi. 

3. bunch w may be associated to an edge of BCj while that edge was already traversed. 

In the first two cases, although w is split and may be part of various boundary-/waveform- 
sections, the invalid marks in the corresponding boundary-/waveform-section trees guarantee that 
these can not cause Type-I/Type-II events. To handle the third case, as soon as we backtrack to 
the time at which boundary cycle BC was split, we are initiating all possible bunches corresponding 
to boundary-sections and waveform-sections which will traverse g in future (as determined by the 
wavefront propagation within REGi). This facilitates in moving a bunch (resp. waveform-section) 
in B forward with another section of wavefront which is also associated with the same boundary- 
section (resp. edge) as B. Hence the third case is not possible. 

The correctness of updating associations, boundary-sections and waveform-sections of a given 
two sections of wavefront relies on the MERGE procedure, whose correctness is proved in the 
Lemma 17.1.11 □ 



Lemma 8.4 The total number of merges during the entire algorithm are 0(m). 
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Proof: The merging is done in a partially traversed junctions or corridors, and there are 0{m) 
corridors or junctions. □ 



Lemma 8.5 The total number of gateways/gateway-edges are 0(m). 

Proof: Since there are a total of 0(m) corridors/j unctions together, the number of possible bound- 
ary cycles are 0(m). ^From Lemma 18. 11 each boundary cycle except the one having t is attached 
with one outgoing gateway. Hence the complexity. □ 

Theorem 8.1 The total number of Type- I /Type- II events due to gateways are 0{m). 

Proof: Let g be an outgoing gateway w.r.t. boundary cycle BC, and let ge be a gateway-edge of g. 
Let REG be the region bounded by BC. A gateway-edge is processed only when all the processing 
is complete for REG. Hence, a side of gateway-edge can cause at most O(l) Type-I/Type-II events. 
Since there are 0(m) gateway-edges (Lemma 18. 5p . the total number of Type-I/Type-II events due 
to gateways are 0(m). Therefore, the total number of Type-I/Type-II events are 0(m). □ 

Lemma 8.6 The amortized time complexity in orienting all the gateways is 0(m(lg m)(lg n)). 

Proof: We can solve point location problem among two polygons, where each is having 0(mi) 
and 0(771-2) edges respectively, in 0(min (mi, 777,2)) time (this approach has been used before in 
[16]). However, we are interested in the case of two polygons bounded with 0(mi) and 0(777-2) 
edges/convex chains with each convex chain having 0(n) edges. The planar point location algorithm 
given in [16] can be extended to yield an algorithm with the time complexity 0(min(mi, 777,2) lg n) 
for this case. For each boundary split, we need to locate either t or a gateway, and, the boundary 
splits occur recursively. Then the recurrence representing this recursion is T{m) = T(mi)-\-T{m2) + 
minjmi, 777-2} lg ?i, where T(m) is the time required to solve the location problem in a polygonal 
region with 0{m) corridor chains and entry/exit boundaries. This results in an overall complexity 
of 0(?77,(lg m)(lg n)). □ 

Lemma 8.7 The total number of Type-Ill events generated in merge procedure are 0(m). 

Proof: We will defer the proof till later where in fact we will show that the total number of Type- 
Ill events in the entire algorithm are bounded. □ 



9 Event Points 

9.1 Type-I and Type-II Events 

The Type-I event determination involves computing the shortest distance for a segment in a 
waveform-section to strike the associated corridor convex chain or enter/exit boundary. The Type- 
II event determination involves finding the shortest distance for the bunch to expand before striking 
either a corridor convex chain or enter /exit boundary in the boundary-section associated with it. 
These shortest distance computations are explained in Section [5j Both the Type-I and Type-II 
events change dB. 

Let J be an untraversed junction i.e., no edge/vertex of it has been traversed. When a section 
of wavefront SW = {Bi, . . . , Bj, . . . , B&} strikes an edge e± = (vi,V2) of J, SW may be split into 
at most two sections: a sequence of wavefront segments SW\ that could possibly strike an edge 
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Figure 25: Split of a section of wavefront 



e2 of J and a sequence of wavefront segments SW2 that could possibly strike an edge e 3 of J 
such that SWi U SW2 = SW. See Fig. l25Ta). The edge e\ in dB is replaced by e<i and e 3 . We 
intend to determine two successive inter-bunch I-curves, l-cwcve(Bj prev , Bj), l-cwve(Bj, Bj gucc ) in 
the wavefront such that l-cmve(Bj prev , Bj) intersects e 2 and l-curve(Bj , Bj gucc ) intersects e 3 . Here, 
Bjprev (resp. Bj succ ) is the bunch that precede (resp. succeed) Bj in the wavefront. This allows 
associating Bi, . . . , Bj-i to ei and -Bj+i, ■ ■ ■ , -Bfc to 63. We do binary search over the inter-bunch 
I-curves of SW to find such j. In this binary search, primitive operation computes the intersection 
of a inter-bunch I-curve with either e2 or e 3 . Since I-curves are higher-order curves, explicitly 
computing their intersections with ei and/or e% is not efficient. Hence we use Ilntersect proce- 
dure described in Section [6j Once we determine such j, we do binary search over the I-curves of 
bunch Bj = (w(vi), . . . , w(v r ), . . . , w(v u )) to compute two I-curves such that I-curve(w(v r _i), w(v r )) 
would intersect e2 and I-curve(u;(w r ), w(v r+ i)) would intersect e 3 . Since intra-bunch I-curves are 
straight-lines, we compute the intersection between an I-curve and an edge. Then using bunch split 
procedure listed in Section 14.14 we split Bj into Bj , Bj such that Bj comprise the wavefront seg- 
ments {w(vi), . . . ,w(v r — 1)}, and B'J comprise the wavefront segments {w(v r ), . . . ,w(v u )}. Then 
we insert w(v r ) to B'- so that w(v r ) is present in both B'- and B" as it can strike either e2 or e%. 
After that, SW\ comprise Bi, . . . , Bj-i, Bj, and SW2 comprise B", Bj+i, . . . , 

Let C be an untraversed corridor i.e., no edge/vertex of it has been traversed. Also, let C 
has CCi,CC2 convex chains, and 61,62 enter /exit bounding edges. When a section of wavefront 
SW = {Bi, . . . , Bj, . . . , Bfr, . . . , Bi} strikes an enter boundary 61 of C, SW may be split into at most 
three sections: a sequence of wavefront segments SW\ that could possibly strike CC\, a sequence of 
wavefront segments SW2 that could possibly strike an edge 62, and a sequence of wavefront segments 
SW 3 that could possibly strike an edge CC 2 of C such that SWi U SW 2 U SW 3 = SW. See Fig. 
[251( b) . The 61 in dB is replaced by CC\, b 2 , CC2. We determine two successive inter-bunch I-curves, 
l-cmve(Bj prev ,Bj), I-curve(Bj, Bj succ ) in SW such that I-curve (Bj prev , Bj) intersects CC\ and I- 
curve(5j, Bj aucc ) intersects 62. Here, Bj prev (resp. Bj aucc ) is the bunch that precede (resp. succeed) 
Bj in the wavefront. Once we determine such j, we do binary search over the I-curves of bunch Bj 
to find w(v r ) G Bj which could possibly intersect both CC\ and 62. As in the case of junction J, we 
split Bj such that w(v r ) £ Bj is in both Bj and B'J. The binary search follows the same procedure as 
listed above (the case of junction J). Similarly, we determine two successive inter-bunch I-curves, 
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l-cmve(Bf !prev ,Bk), l-cuvve(B k , B ksucc ) in SW such that l-cuive(B kprev , B k ) intersects b 2 and I- 
cmve(B k ,B ksucc ) intersects CC 2 and split B k into B' k and B' k . Here, B kprev (resp. B ksucc ) is the 
bunch that precede (resp. succeed) B k in the wavefront. Then SW± comprise {Bi, . . . , Bj^%, Bj}, 
SW2 comprise {Bp Bj+%, . . . , -Bfc-i, B' k }, whereas SW3 comprise {B' k , B k+ \, . . . , Bi}. Due to the 
non-crossing nature of I-curves (Lemma 12. ip . point of tangency corresponding to Type-Ill event on 
CC\ due to SW\ is caused by the wavefront segment w{v r ). The distance to point of tangency p 
on CC\ along v T p from the periphery of w(v r ) is pushed to the event heap. Another Type-Ill event 
due to the interaction of CC 2 and SW3 is also pushed to the min-heap. As a whole two Type-Ill 
event points, and a shortest distance between SW% and b 2 are the event points determined during 
this procedure. Note that no other interactions between SW\ and CC\ or SW3 and CC 2 worth 
further consideration, as the destination t is in its own corridor. 

Let a bounding edge of J (resp. C) was traversed, and the wavefront struck the other bounding 
edge of J (resp. C). This causes boundary split. The merging procedure is explained in section [71 

Consider a junction J = (ei, 62,63) that is traversed. Suppose the edge e\ is traversed whereas 
the edges e2 and ez are not traversed when W(d). Also, suppose the bounding edge e2 is struck 
when W(d') This causes boundary split. While updating the relevant associations of sections of 
wavefront, the merging procedure described in Section [7J may in turn cause Type-Ill events. The 
same is true when a boundary splits due to merger in corridors. 

A section of wavefront SW striking a junction (resp. corridor) boundary vertex is handled 
similar to the above cases, except that both the edges adjacent to that vertex are considered as 
struck by the wavefront. 

Lemma 9.1 The total number of bunch/waveform-section splits during the entire algorithm are of 
0(m). 

Proof: The waveform-sections are associated with the edges which are not struck as of now. The 
split occurs when a waveform-section strikes the boundary of a junction/corridor and there are 
0(m) such boundary edges in all the junctions/corridors together. □ 

Lemma 9.2 The complexity of bunch/waveform-section splits during the entire algorithm are 
0(m(\g m)(lg n)). 

Proof: Each split takes O(lgn) time as it involves binary search over the I-curves. The appropriate 
boundary- and waveform-section updates take 0((lgm)(lgn)) time (Lemmas I4.5(l4~6|) . Combining 
this with Lemma l9.lt leads to the proof. □ 

Lemma 9.3 The total number of Type-I and Type-II events together are 0{m). 
Proof: The Type-I or Type-II events occur due to the following: 

(1) A section of wavefront strikes either a corridor convex chain or a corridor enter /exit boundary. 

(2) A section of wavefront striking gateway-edges in merge procedure. 

Since there 0(m) corridors, there are 0(m) Type-I/Type-II events due to (1). ^From Lemma |8.1[ 
the number of Type-I/Type-II events due to (2) are bounded with 0(m). □ 

Lemma 9.4 The total time in determining and handling all Type-I/Type-II events is 0(m(lgm)(lgn)). 
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Proof: ^From Lemma 15.21 determining Type-I event is of complexity 0((lgm)(lgn)). /,From 
Lemma l5.1j, determining Type-II event is of complexity 0((lgm)(lgn). 

Both of these event types involve updating WST and/or BST, either by splitting or by invok- 
ing the MERGE procedure. ^From Lemma 14.51 updating WST takes 0((lgm)(lgn)) amortized 
time. ^From Lemma 14.61 updating BST takes 0((lg m)(lg n)) amortized time. /.From Lemma 18.61 
amortized time in orienting all the gateways is 0(m(lgm)(lgn)). ^From Theorem 111.31 amortized 
time involved in the merge procedure is 0((lg m)(lg n)). 

Type-Ill event determination and handling costs are considered in 111.11 

□ 



9.2 Type-Ill Event 

This event occurs when a wavefront segment w(v r ) in a section of wavefront SW strikes a corridor 
convex chain CC at a point of tangency p along the tangent v r p from v r to CC. The event 
determination procedure is listed as part of Type-I and Type-II event handling. See Subsection 
19. 1L The event possibly causes either the initiation of a new bunch from p, or modifying an 
existing bunch initiated from a vertex of the convex chain CC. The cases in which a bunch is 
initiated/modified are explained in Section [4.11 

Lemma 9.5 The total number of Type-Ill events are O(m). 

Proof: Since there are 0{m) corridor convex chains, number of BHTs initiated (in Case (1) of 
Subsection I4.ip are 0(m). Let C = {v\,V2, • • • v n >} be a convex chain from which a bunch is ini- 
tiated. New bunches from the same convex chain may be introduced during the merge procedure 
(corresponding to Case (4) of Subsectior i4.ip . We note that at the introduction of a new bunch hull 
tree an old bunch of segments, call it B\ is removed from the wavefront. Thus at most one bunch 
B{v z ,v n i exits the bounding edge of the corridor defining the chain C. Let B p be the bunch that 
led to the generation of the bunch B\. The removal of bunches can thus be charged to B p . We thus 
need to bound the number of bunches that a particluar bunch, say B p can generate. If a bunch B p 
generates more than two bunches, then it does so in two different corridors. In this case, the bunch 
B p is split into two at a vertex vj in junction J. The shortest path to this junction vertex is thus 
determined due to a segment in B p (due to the non-crossing property of bunches). Thus the split 
of bunches can be charged to junction vertices. These are 0{m) in number. And each bunch that 
is generated and removed, as B\ is above, can be charged to a bunch or split portion of a bunch. 
This gives the desired bound of 0(m). □ 



9.3 Type-IV Event 

The intersection of inter-bunch I-curves within a waveform-section are captured with this event. 
Using these points of intersection, we can detect when two non-adjacent bunches within a waveform- 
section meet. Since the intra-bunch I-curves are diverging, only the intersection of inter-bunch 
I-curves and I-curves from two different bunches are considered. 

Lemma 9.6 Let d! be the shortest distance between a waveform- section WS and its associated ver- 
tex/edge e, which corresponds to an event evt. Even though some of the I-curves in that waveform- 
section intersect among themselves before evt occurs, WS does strike e for the first time when the 
evt occurs (provided that the association of e does not change due to mergers). 

Proof: 

See Fig. [26j Suppose the inter-bunch I-curves in WS, l-cwve(w(ui),w(u2)) and l-cuive(w(vi),w(v2)), 
intersect at point p, whereas the Euclidean distance from p to WS is d", for d" < d' . Consider the 
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(a) A view of I-curve intersection in a W ST 



(b) Event point distance before and after I-curve in- 
tersection 



Figure 26: I-curve intersection and the computed shortest distance 



wavefront just after the evt occurred i.e., after it traverses an Euclidean distance e (for a positive 
constant e). At this wavefront progression, let WS' be the waveform-section comprising bunches 
in WS sans the the bunches removed due to the inter-bunch I-curve intersection. Let w'{u\) (resp. 
w'(u2), w'{v i), w'{v2)) be the wavefront segment in WS' corresponding to the wavefront segment 
w{u\) (resp. w(u2),w(vi),w(v2)) in WS. Consider any point q s.t. q is external to both WS and 
WS', and, q is closer to w(ui) than w{v2) (symmertic cases can be argued similarly). Let x(y) 
be the point of intersection of line segment qui (resp. qv2) with w(ui) (resp. w(v2))- Given that 
\\qx\\ < \\qy\\, it is trivial to note that (\\qx\\ — (d" + e)) < (\\qy\\ — (d" + e)). In other words, it 
is guaranteed that the point q is closer to w'(ui) than w'{v2) whenever q is closer to w(ui) than 
w(g 2 ). □ 

However, for the sake of utilizing I-curves for further computations, and to avoid overlap of 
wavefront segments within the wavefront, we detect their intersection. There are three types of 
I-curve intersections possible in a waveform-section: 

(1) Intersection of adjacent inter-bunch I-curves 

(2) Intersection of non-adjacent inter-bunch I-curves 

(3) Non-adjacent I-curve intersection causing partial elimination of bunches 

The dirty bridges discussed in Section [10] takes care of Case ([3]). Since explicitly computing the 
I-curve intersections occurring within a WS is compute intensive, both the I-curve intersections 
mentioned in Case (pQ) and Case ([2]) are detected by finding the shortest distance between the 
sibling hulls stored at the internal nodes of WST. Let UHi,UH r be the hulls at sibling nodes vi 
and v r of a WST respectively. Let b\ m (resp. b 1 ™) be the rightmost (resp. leftmost) bunch in the 
bunches stored at the leaves of vi (resp. v r ). Also, let UH[ m (resp. UH l r m ) be the hull of b\ m 
(resp. &[ m ). Since two adjacent bunches always intersect along an I-curve, we avoid detecting event 
corresponding to the intersection of adjacent bunches. Hence we compute the shortest distances' 
d', d" between the hulls UH t - UH[ m , UH r and UH h UH r - UH l r m respectively. And, we push the 

Type-IV event point with distance min(y, 4j-) to the event heap. 

Lemma 9.7 Two I-curves in a WST intersect if and only if Type-IV event occurs. 

Proof: Let bi,bj be two bunches occurring in the left-to-right ordering of bunches in a waveform- 
section. Let bi gucc is the successor of 6« in the wavefront; also, let bj d is the predecessor of bj in 
the wavefront. Suppose the I-curve bi succ ) and I-curve(6j pred , bj) intersect at point p. ^From the 
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(a) Adjacent I-curve intersection 



(b) Hull intersection 



Figure 27: Type-IV Events 

I-curve definition, p € (6j PI bj) when the bunch bi or bj strikes p. Consider the least common ances- 
tor node v of nodes that store bi and bj. The left and right children of v, say v\ and v r , implicitly 
store upper hulls UH[, UH r such that one contains 6, and the other contains bj. Let \JH'[ m (resp. 
XJH 1 ™)) be the hull of the right-most (resp. left-most) bunch in the bunches stored at the leaves 
of vi (resp. v r ). It is immediate to see that the bunches bi,bj intersect, whenever either the hulls 
UHi - UH[ m , UH r or the hulls UH h UH r - UH l r m intersect. This is captured as a Type-IV event. □ 

As shown in the following Lemma, a Type-IV event causes the disappearance of one or more 
bunches, in turn, the waveform-section is dynamically updated. Using the Ilntersect procedure 
(detailed in Section [6]) with RV of a section of wavefront and an I-curve as parameters, we determine 
new associations and update the shortest distance in these new associations. As part of this update, 
we compute the new inter-bunch I-curves among the bunches which became adjacent and their 
intersection points. The event point min-heap is updated accordingly. The following Lemma is 
useful in Analysis. 

Lemma 9.8 Whenever a Type-IV event occurs, there exists at least one bunch which is not required 
to be progressed further. 

Proof: Consider two adjacent inter-bunch I-curves, I-curve(6j_i, bi) and I-curve(6j, h + i). Suppose 
the Type-IV event occurred due to the intersection of these two I-curves (Fig 127(a)). Let r be the 
point of intersection. Also, let REG be the untraversed region bounded by these two I-curves and 
the bunch bi. Suppose the destination t is located in REG. Then due to the no n-cross ing property 
of shortest paths (Lemma l2.ip . the bunch bi is not required to be progressed in REG. Further the 
bunches other than bi-\,bi and bi + \ need not be propagated further. Suppose the destination t is 
located in the region REG. Again due to the non-crossing property of shortest paths (Lemma l2.ip . 
it is not required to progress the bunch bi (in REG). 

Suppose the Type-IV event occurred due to an event point that corresponds to the distance 
between two hulls and does not correspond to the intersection of adjacent I-curves. When this 
event occurs, let CH'^CH^ be the hulls at nodes I and r of WST, and, p be the point at which 
the hulls CHi,CH' r intersect. Let B[ (resp. B' r ) be the set of bunches at all the leaf nodes of the 
subtree formed with the descendents of node I (resp. r). We define the bunches bi and bj herewith: 

• Suppose the point p £ b, for any bunch b S B[ i.e., the point p is located on a bridge of CH[. 
Let B" C B', be the set of bunches, where a bunch b € B" if and only if there exists a line 
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W(Vj) V J 



Figure 28: Dirty Bridges in WST 

segment LS from p to a point located on b s.t. the interior of LS does not intersect with any 
bunch. Among all the bunches in B", we let the bunch stored at the rightmost leaf node of 
WST be termed as 6j. Otherwise, suppose the point p is located on a bunch b S B[. Then 
the bunch b is termed as b{. 

• Suppose the point p £ b, for any bunch b € B' r i.e., the point p is located on a bridge of CH' r . 
Then the set B" C B' r is defined similar to the above case. Among all the bunches in B", we 
let the bunch stored at the leftmost leaf node of WST be termed as bj. Otherwise, suppose 
the point p is located on a bunch b G B' r . Then the bunch b is termed as bj. 

Let B be the set of bunches strictly between 6j and bj in the waveform-section. B is non-empty 
since the event corresponds to an intersection of non-adjacent I-curves. Let the set S comprise 
untraversed regions s.t. a region R £ S if and only if the closed set representing the region R 
intersects with at least one bunch in B. Also, let the region REG be the union of the regions in S. 
Note that the region REG is not necessarily connected. 

Suppose the destination t is located in REG. Let p' be any point chosen on any bunch b in B. 
For any point p" in REG, the line segment p'p" intersects the boundary of CH', U CH' r . In turn, 
p'p" intersects a bunch contained in CH! U CH' r . Then from the non-crossing property of shortest 
paths (Lemma 12. ip . bunch b is not required to be progressed further. In other words, no bunch 
belonging to B need to be progressed in REG. 

Suppose the destination t is located in REG. Then by defining bunches b{,bj symmetric to the 
above, the proof would be similar to the case mentioned in the above paragraph. □ 

As required in the above Lemma, when a Type-IV event occurs we determine the region in 
which the destination point t resides. This is done using the procedure from Vaidya [16j . 

10 Dirty Bridges in WST 

Consider the maintenance of bridges at internal nodes as the wavefront progresses. For pj € 
w(vj),pk € w(vk), let br p = (pj,Pk) be a bridge at node p of a WST when W(d), and Pj,Pk are 
respectively the points of tangencies of w(vj) and w(vk). Let p'j be a point at Euclidean distance 
(\vjPj\ + d! — d) from Vj along the vector pjv). Also, let p' k be a point at Euclidean distance 
{vkPk + df — d) from Vk along the vector pk^l- As the wavefront expands, the bridge br p moves as 
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Figure 29: Dirty bridges versus the correctness of shortest distance computations 



fast as any point on a bunch b r , for i < r < k, expands (Property I4.3p . Therefore, the traversal 
of pj (or pk) could happen due to some bunch b r , where r < i or r > k i.e., b r would be located 
in a subtree other than the one at p. Let an endpoint of br' p , say p'j, does not belong to W(d') for 
d! > d. Since the Property 14.31 is no more applicable, br' p is not a bridge at node p at W(d'). As 
shown in Fig. [28], due to the definition of a bridge in WST, br' p ' is the bridge at node p when W(d'). 
However, it is inefficient and impractical to update the bridges infinitely often at nodes such as p. 
As long as at least one wavefront segment in each of the bunches bj and bk (on which p'- and p' k 
incident) belongs to W(d') i.e., as long as bj and bk are at the leaves of WST corresponding to 
W(d'), we continue to have br p as a dirty bridge at node p of WST. 

Definition 10.1 A bridge is considered as a dirty bridge whenever it is defined at a WST node, 
but one or both of its endpoints do not belong to the shortest path wavefront. If the bridge is dirty 
at a node p of WST, then the corresponding upper hull at p is termed as a dirty upper hull. 

When bridge br p is determined, say at W(d), at a node p of WST, let p'j 6 w(vj) and pk G w(vk) 
be its endpoints such that p'j is not a point of tangency of w(vj). Let C be a circle with the same 
center and radius as w(vj). To facilitate in reconstructing bridge at W(d'), by utilizing property 
14.31 we find a point of tangency pj on C (closest along C to p'j) so that PjPk is a tangent to both 
w(vj) and w(vk). Again, so formed bridge br p = (j)j,Pk) is a dirty bridge. 

Lemma 10.1 The dirty bridges in a WST do not affect the correctness of shortest distance com- 
putations between WST and its associated corridor convex chain or enter/exit boundary C . 

Proof: Consider a node p of W ST with the bunches at its descendant leaf nodes as hi, . . . ,bj, 
. . . , bk, ■ ■ • , b w -\,b w . See Fig. [29j Let B = Uig{i to} Let a dirty bridge br p at node p of WST 
join two points p'j € bj and p' k € bk- Suppose the point pj was traversed by a bunch b u , where 
u $l {i, i + 1, . . . ,w}, when W(d') . Let ch' p be the dirty upper hull at node p. Consider bridge 
br' q at node p computed at this stage of wavefront progression. Since we know that at least two 
bunches bj and bk exist among bunches in B, the bridge br' q is guaranteed to exist. Let the bridge 
br' q joins two points p' r and p[, where the point p' r is located on the wavefront segment w(v r ) £ b r 
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whose center is v r and the point p\ is located on the wavefront segment w(vi) G b\ whose center is v\. 
Note that w(v r ) (resp. b r ,p[, w(vi), bi) is not necessarily distinct from w(vj) (resp. bj,p' k ,w(vk),bk)- 

Consider a wavefront progression W(d) at which the bridge at node p was not dirty. Let br p 
be the bridge at node p for W(d) with Pj,Pk as its endpoints. Let the corresponding upper hull 
be ch p . Choose a point p r £ w(v r ) at that stage of wavefront expansion so that the line segment 
v r p r extended yields p' r . Then the Euclidean length of line segment p r p' r be d! — d. /,From the 
property 14,31 for an expansion of d' — d along v r p' r , the Euclidean length of segments PjPj, pw[ 
would be d! — d too. This with the fact that the bunches remaining in B after a Type-IV event are 
in the same order along B as they were before assures the following: p' r 6 ch' p and p\ € ch' p . The 
convexity of ch' p guarantees that the line segment br' q (joining points p' r ,p'i) belongs to ch' p . Hence 
the shortest distance from C to br' p is less than (or equal to) the shortest distance from C to br' q . 
Similar arguments can be given when both the endpoints of a dirty bridge were traversed. 

The correctness of the shortest distance computation procedure relies on the property shown 
above for a dirty bridge, when the shortest distance from C is to that dirty bridge. Suppose the 
shortest distance between br' p and C is less than the shortest distance to either of the bridges at p's 
immediate descendants. Since the shortest distance from C to bridge br' p is less than (or equal to) 
the shortest distance from C to bridge br' q , the bridge br' p will be struck whenever the wavefront 
can strike the bridge br' q . Therefore, the correctness of shortest distance computation procedure is 
not affected. □ 



Lemma 10.2 The bridges constructed over invalid segments of a BHT do not affect the correctness 
of shortest distance computations. 

Proof: The argument is similar to the last paragraph of Lemma 110.11 □ 



Lemma 10.3 The total number of Type-IV events are 0{m). 

Proof: The Type-IV events are caused by those bunches which passed over the exit boundary of 
the corridor in which they were initiated. Since there are 0(m) such bunches (Lemma 19. ip . there 
can be 0{m) Type-IV events. □ 

Lemma 10.4 The total cost in determining and handling all the Type-IV events is 0(m(lgm)(lgn)). 

Proof: Initially, in all WSTs together, there are 0{m) bunches. The insertion or deletion of a 
bunch B to a WST affect only the Type-IV events along the path from the leaf node at which 
B is inserted to the root. Since the depth of a WST is 0(m), computing the shortest distance 
between two upper hulls in WST is O(lgre), the amortized time to compute Type-IV events and 
update the min-heap takes 0((lg m) (lg n)) time. Only the bunches that traversed the exit boundary 
of the corridor in which they are initiated would require these updates, which are 0(m) in number. 

As mentioned in Lemma 19.81 each Type-IV causes at least one bunch to be removed from the 
given WST. ^From Lemma 1431 updating WST corresponding to this deletion takes 0((lg m)(lg n)) 
amortized time. Also, from LemmaE21 updating shortest distances in the min-heap takes 0((lg m)(lg n) 
amortized time. Since there are at most 0{m) Type-IV events, total computation involved in up- 
dating WST takes 0(m(lgm)(lgn)). 

The amortized time complexity involved in locating the destination t using the procedure from 
|16j takes 0((lgm)(lgn)) time (Analysis is similar to Lemma 18.61 ). The updates to WST and 
the shortest distance computations due to a Type-IV event take 0((lgm)(lgn)) time. The dele- 
tion/updating of events in the event queue take O(lgre) amortized complexity. □ 
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11 More Analysis 

Theorem 11.1 The total number of event points are 0(m). 
Proof: ^From Lemmas 19.31 19.51 and 110.31 



□ 



Theorem 11.2 The number of waveform- and boundary-sections created/updated are 0{m). 

Proof: Since there are at most O(m) bunches at any point of execution of the algorithm (Lemma 
14.41) and a bunch is part of at most one boundary-section, initially there are at most 0(m) boundary- 
sections. Since there are 0(m) corridor convex chains/exit boundaries and each can have at most 
one waveform-section associated to it, initially there can be at most 0(m) waveform-sections. We 
need to update these initial boundary- and waveform-sections due to events, which are upper 
bounded by 0(m) (Theorem lll.ip . Hence the complexity. □ 

Theorem 11.3 The total processing involved in the merge procedure, excluding the computation 
involved in finding Type-Ill events, during the entire algorithm is of 0(m(lg m)(lg n)) time com- 
plexity. 

Proof: We need to invoke MERGE procedure constant number of times at most per each merge 
operation. The MERGE procedure complexity relies on ASSOCATOB procedure. Primarily, the 
complexity of the ASSOCATOB procedure relies on two components: how many times we call the 
binary search procedure within it, and on the total number of WST, BST updates. 

In the procedure, the binary search over the I-curves of B determines the next I-curve which 
does not intersect the edge c prev . The RV of a bunch may be before a merge, or it can become 
after calling the ASSOCATOB procedure. However, the latter is possible only when the merge 
occurs and the number of merges are bounded by 0{m) from Lemma 18.41 ^From Theorem lll.2[ 
there can be at most 0(m) waveform-sections during the entire algorithm, and hence there can be 
at most 0(m) bunches with RVs. The RV of a bunch can change from to non-0 because of 
splits, however, the splits are bounded by 0(m) (from Lemma \9.1\i : hence there are 0(m) bunches 
with non-0 RVs over the course of the algorithm. Since each binary search over the bunches finds 
a bunch whose RV is non-^, we charge the O(lgm) binary search complexity to bunches with 
non-0 RVs, leading to O(mlgm) complexity. The Ilntersect procedure takes O(lgm) time per 
invocation. Also, determining the proximity of a corridor convex chain or enter /exit boundary c s 
(line 4 of ASSOCATOB procedure) w.r.t. a bunch or WST in A versus a bunch or WST in B 
takes 0((lgm)(lgn)). Hence, the overall time is 0(m(lg m)(lg n)). 

Since there are 0{m) waveform and boundary-sections possible in the entire algorithm (Theorem 
111. 2D . and, each WST/BST update takes 0((lgm)(lgn)) (using [11]), we spend 0(m(lgm)(lgrz)) 
to dynamically maintain waveform- and boundary-sections during all the invocations of ASSOCA- 
TOB procedure together. 

Also, assigning all the corridor convex chains or enter/exit boundaries prior to the intersection of 
the first inter-bunch I-curve(&i, 62) of B with the RV(A), and, adjusting the associations of already 
existing corridor convex chains/exit boundaries in RV(b\) takes 0((lgm)(lgn)) time. Hence we 
spend 0(m(lgm)(lgn)) time in ASSOCATOB procedure during the entire algorithm. Therefore, 
MERGE procedure, and, all merge operations together are of 0(m(lgm)(lgn)) time complexity. □ 

Lemma 11.1 Building and maintaining all BHTs during the entire algorithm takes 0{m\gn + n). 
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Proof: The BHT maintenance is divided into four cases (See Section l4.1i ). In Case (1), we create 
nodes corresponding to the bunch vertices all at once in 0(n) time, even though some of them are 
invalid by the time we construct BHT. In Cases (2) and (3), we do nothing. In Case (4), splitting 
and balancing T new takes O(lgn) time. Since there are 0(m) waveform- and boundary-sections 
(Theorem II 1.2|) together, Case (4) takes O(mlgn) time in the worst-case. □ 

Theorem 11.4 Any event point determination with the associated updates, excluding the initiation 
and maintenance of BHT s, can be accomplished in 0((lgm)(lgn)) amortized time. 

Proof: Immediate from Lemmas 19. 4[ 111.11 and 110.41 □ 

Theorem 11.5 The shortest distance from s to t can be computed in 0(T + m(lgm)(lgn)) time 
using 0(n) space. 

Proof: The triangulation of polygonal region takes 0(n + mlg 1+<E m), represented as 0(T). Given 
a triangulation, finding useful corridors and junctions takes O(mlgn). ^From theorems 111.11 Hi. 41 
all the event points determination and associated updates, excluding the updates to BHTs, can 
be done in 0(m(lgm)(lgn)) time. ^From Lemma 1 11. 11 initiating and maintaining all BHTs take 
0{m\gn + n). Hence the algorithm is of 0(T + m(lgm)(lgn)) time complexity. 

All the four data structures BHT, BST, WST and Event Heap require at most 0(n) space at 
any instance during the entire algorithm. Hence, the algorithm is of 0(n) space complexity. □ 



12 Conclusions 

We have described an algorithm for finding the Euclidean shortest path in polygonal domain with 
0(T + m(lgm)(lgn)) time complexity using 0(n) space. It would be of interest to investigate for 
a solution with 0(n + mlgm) time and 0(n) space. Also, exploring the applicability of the above 
technique to weighted geodesic shortest path computation, and determining approximate Euclidean 
shortest paths in polygonal domains is of interest. 
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